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1.0  ABSTRACT 


The  investigation  of  an  aeroelastically  induced  constrained  warping 
Phenomenon  for  a  comoosite,  suoermaneuverable  tyne  aircraft  wing 
has  continued  in  this  second  year  of  the  study.  The  first  year 
investigation  was  concentrated  mainly  on  the  static  phenomena  and 
the  search  for  closed  form  solutions  for  free  vibration  of  air¬ 
craft  wings  having  constrained  warping  in  the  presence  of  elastic 
counting.  The  wing  is  analytically  modelled  as  a  straight  flat  lam¬ 
inated  Plate.  Various  forms  of  highly  simnlified  aerodynamic  loads 
are  emoloyed  in  the  analysis.  The  free  vibrations  and  stability 
aspects  of  this  phenomenon  are  examined  to  obtain  some  physical  in¬ 
sights  and  to  determine  its  importance  and/or  design  imnlications . 
Analytical  tools  employed  include  an  affine  transformation  concent 
which  was  formulated  previously  (by  the  present  principal  investi¬ 
gator)  as  well  as  a  non-dimensionalization  scheme.  With  the  helm 
of  these  tools,  an  evolution  of  effective  warping  narameters  with 
which  to  study  this  phenomenon  was  carried  out.  The  virtual  work 
theorem  and  variational  principles  were  used  to  derive  the  equa¬ 
tions  of  motion  based  on  the  assumed  wing  disnlacements .  The  exis¬ 
tence  of  closed- form  free  vibrations  solutions  for  composite  wincrs 
with  elastic  coupling  and  constraint  of  warping  was  established. 

A  revelation  from  these  closed-form  solutions  is  that,  elastic 
couplinq  lowers  the  couoled  frequencies  (in  fact  a  significant  a- 
mount  of  coupling  could  reduce  the  first  freauency  to  almost  zero, 
hue  discovery  of  the  closed- form  solutions  for  the  free  vibration 
which  seems  to  mark  the  first  time  such  solutions  were  ever  ob¬ 
tained,  not  only  led  to  answers  to  a  number  of  previously  unan¬ 
swered  nuestions  but  also  raised  new  unanswered  questions  such  as 
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"Does  the  aeroelastic  divergence  of  flutter  problem  of  such  winqs 
have  any  closed-form  solutions". 

A  main  goal  of  the  second  year  investigation  is  therefore  essen¬ 
tially  to  find  the  answer  to  the  above  mentioned  question.  Accor¬ 
dingly  in  the  first  quarter,  efforts  were  concentrated  on  deter¬ 
mining  the  possibility  of  finding  closed-form  solutions  to  the  di¬ 
vergence  problem.  The  investigation  led  to  two  possible  methods 
of  obtainino  such  solutions.  These  methods  are  (a)  the  "elimina¬ 
tion"  approach  (which  was  used  for  the  free  vibration)  and  (b)  the 
method  of  Laolace  transformation.  Although  it's  already  known 
that  (a)  works,  (b)  was  implemented  in  formulating  the  analytical 
expressions  for  the  closed-form  solution  of  the  divergence  prob¬ 
lem  due  to  an  anticipated  relative  ease.  These  newly  derived 
closed- form  expressions  are  shown  in  this  report.  The  numerical 
a?proach  used  for  extracting  the  eigenvalues  from  the  closed-form 
expressions  for  the  free  vibration  was  studied  carefully  to  de¬ 
termine  how  it  might  be  used  to  extract  the  eigenvalues  (and 
eventually  the  divergence  flight  velocities)  from  the  closed  form 
expressions  derived  for  the  divergence  phenomenon.  The  results 
extracted  showed  agreement  with  existing  results  and  seem  to  es¬ 
tablish  new  trends. 
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2.0  NOMENCLATURE 


=  chordwise  integrals 

=  affine  space  half-chord  and  chord,  respectively 
=  elastic  constants 
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parameter  that  measures  the  location  of  the 
reference  axis  relative  to  mid-chord 

bending  and  torsional  stiffness,  respectively 

wing  box  depth 

affine  space  bending  and  torsional  displacement, 
respectively 

elastic  coupling  and  warping  stiffness,  respec¬ 
tively 

elemental  stiffness  parameter 
Strouhal  number 

affine  space  running  aerodynamic  lift  and  mo¬ 
ments,  respectively 

affine  space  half-span  for  the  wing 

affine  space  mass  per  unit  span 

differential  aerodynamic  pressure  distributions 
in  physical  and  affine  space,  respectively 

time 

virtual  work  expressions  in  physical  and  affine 
space,  respectively 

flutter  speed 

physical  and  affine  space  coordinates,  respec¬ 
tively 

displacement  shape  functions 

generic  nondimensionalized  stiffness  parameters 
affine  mass  ratio  parameter 

Poisson  ratios  and  generalized  Poisson's  ratio, 
respectively 

divergence  parameter 
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re- 


affine  space  material  and  air  density, 
spectively 

twisting  displacements 
displacement 
vibration  frequency 
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3 . °  INTRODUCTION 


In  the  late  1960's  designers  began  to  investigate  the  possibility 
of  exploiting  the  directional  properties  of  composite  materials 
to  improve  aeroelastic  characteristics  of  lifting  surfaces.  For 
example,  superfighter  designers  may  aeroelastically  tailor  a  wing 
so  that  it  deforms  to  the  optimum  camber  under  maneuvering  loads. 

For  instance,  a  wing  can  be  tailed  so  that  its  leading  edge  will 
twist  downward  under  the  stress  of  a  tight  turn,  thereby  decreas¬ 
ing  the  wing's  angle  of  attack  and  hence  reduce  drag. 

Exploitation  of  the  directional  properties  of  composite  materials 
to  solve  aeroelastic  instability  problems  (known  as  aeroelastic 
tailoring)  may  enhance  the  performance  of  future  high  performance 
military  aircraft,  particularly  the  super -maneuverable  concepts. 

This  stems  from  the  need  for  future  military  weapon  systems  to  ex¬ 
hibit  high  performance  and  minimum  vulnerability  (e.g.;  minimum 
radar  cross-sectional  area) .  These  requirements  may  lead  to  unor¬ 
thodox  aircraft  configurations  which  in  turn  result  in  unorthodox 
aeroelastic  instability  problems  (e.g.  the  X-29  primary  mode  of 
instability,  called  the  'body  freedom'  flutter). 

Perhaps  the  aeroelastic  tailoring  concept  would  have  been  dis¬ 
covered  much  earilier  if  the  physics  of  anisotropic  aeroelasi icity 
were  more  apparent.  This  inherent  physical  intractability  is  large¬ 
ly  due  to  the  existence  of  numerous  variables,  e.g.;  flight  param¬ 
eters,  several  composite  directional  properties,  fiber  orientation 
angles,  etc.,  which  are  not  even  truly  independent,  in  the  aniso- 
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tropic  system.  This  state  of  affairs  is  therefore  analogous  to 
that  which  existed  in  the  basic  rigid  body  aerodynamics  before 
the  advent  of  the  similarity  rule  theory.  This  theory  clearly 
revealed  that,  the  utilization  of  non-dimensionalized  groups  of 
variables,  e.g.;  Reynolds,  Mach,  Strouhal,  Froude  numbers  etc., 
provide  significantly  superior  physical  insight  to  the  problem 
than  the  individual  physical  variables.  The  new  methodology  that 
is  being  used  as  the  basic  tool  in  this  research  program  is  basi¬ 
cally  the  aeroelastic  equivalent  of  the  aerodynamic  similarity 
rule.  The  expected  superiority  of  this  new  approach  over  the 
state-of-the-art  (SOA)  counterpart,  which  utilizes  individual 
physical  variables,  especially  in  terms  of  physical  insights,  has 
been  demonstrated  in  references  1-10.  For  example,  a  high-aspect- 
ratio  composite  wing  could  behave  aeroelastically  like  a  low  as¬ 
pect  ratio  wing  and  vice-versa.  Similarity  parameters  can  expose 
conditions  for  which  this  might  happen.  This  is  significant,  (for 
instance)  in  the  light  of  the  important  role  played  by  the  wing 
aspect  ratio  in  the  aerodynamics  approximations  for  an  aeroelastic 
analysis.  This  result  may  therefore  be  suggesting  some  new  form 
of  coupling  between  the  elastic  and  aerodynamic  equations  in  com¬ 
posite  wing  aeroelas ticity .  A  fundamental  aspect  of  this  new 
methodology  has  been  used  in  studies  at  Purdue  University- and 
MlT12. 

In  this  research,  an  investigation  of  a  wing's  spanwise  sectional 
distortion  (warping)  resulting  from  aerodynamic  forces  and  its  in¬ 
fluence  on  the  wing's  free  vibrations,  as  well  as  (aeroelastic) 
stability,  is  being  carried  out.  In  particular,  St.  Venant's  tor- 
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sional./ twisting  theory,  which  is  currently  widely  used  for  esti¬ 
mating  the  wing's  twisting  displacements  has  been  examined  with 
the  help  of  the  new  methodology,  to  determine  its  limitations, 
when  applied  to  wings  fabricated  of  composite  materials.  The 
relevance  and  significance  of  this  study  for  the  newly  emerging 
supermaneuverable  type  aircraft  may  be  seen  in  the  light  of  the 
fact  that  (a)  supermaneuverability  is  characterized  by  high  angle 
of  attack  which  implies  high  twisting  aerodynamic  forces,  and 
(b)  most  aircraft  designers  believe  future  aircraft  will  be  fabri¬ 
cated  of  40-70%  composite  materials. 

When  the  St.  Venant ' s  torsional  theory  is  used  to  estimate  a  wing's 
twisting  displacement  and/or  f orces ,  the  fact  that  the  wing's  root 
section's  distortion  is  relatively  small  compared  to  that  of  other 
sections  is  ignored.  However,  previous  investigators  have  deter¬ 
mined  that  such  an  unrealistic  assumption  may  lead  to  only  little 
errors  if  the  aspect  ratio  of  the  wing  is  very  high. 

The  research  has  shown  during  the  first  phase  that  the  conclusions 
reached  by  previous  investigators  are  basically  true  for  wings  fab¬ 
ricated  of  metals  or  isotropic  materials.  A  set  of  new  theories  is 
therefore  being  postulated  in  this  research  effort  for  accurately 
estimating  the  twisting  displacements,  vibrational  frequencies  and 
instability  boundaries  for  wings  fabricated  of  composite  materials. 


7 


REFERENCES 


1.  Oyibo,  G.A.,  "A  Unified  Flutter  Analysis  for  Composite  Aircraft 
Wings",  AIAA  paper  No.  84-0906,  25th  AIAA  SDM  Conference,  Palm 
Springs,  CA,  May  1984. 

2.  Oyibo,  G.A.,  "Reply  by  the  Author  to'p.A.A.  Laura",  AIAA 
Journal,  Vol.  22,  April  1984. 

3.  Oyib->,  G.A.,  "A  Unified  Flutter  Theory  for  Composite  Aircraft 
Wings"  Invited  Lecture,  26th  Israel  Annual  Conference  on  Avia¬ 
tion  and  Astronautics,  Technion,  Haifa,  Israel,  February  1984. 

4.  Oyibo,  G.A.,  and  Brunelle,  E.J.,  "Vibrations  of  Circular  Ortho¬ 
tropic  Plates  in  Affine  Space",  AIAA  Journal,  Vol.  22,  October 
1984  . 

5.  Oyibo,  G.A.,  "Generic  Approach  to  Determine  Optimum  Aeroelastic 
Characteristics  for  Composite  Forwardswept  Wing  Aircraft" ,  AIAA 
Journal,  Vol.  22,  No.  1,  January  1984,  pp.  117-123. 

6.  Oyibo,  G.A.,  "Unified  Aeroelastic  Flutter  Theory  for  Very  Low 
Aspect  Ratio  Panels",  AIAA  Journal,  Vol.  21,  No.  11,  November 
1983,  pp.  1581-1587. 

7.  Brunelle,  E.J.,  and  Oyibo,  G.A.,  "Generic  Buckling  Curves  for 
Specially  Orthotropic  Rectangular  Plates",  AIAA  Journal,  Vol. 
21,  August  1983,  pp.  1150-1587. 

8.  Oyibo,  G.A.,  "Unified  Panel  Flutter  Theory  with  Viscous  Damp¬ 
ing  Effects",  AIAA  Journal,  Vol.  21,  May  1983,  pp.  767-773. 

9.  Oyibo,  G.A.,  "Flutter  of  Orthotropic  Panels  in  Supersonic  Flow 
Using  Affine  Transformations",  AIAA  Journal,  Vol.  21,  February 
1983,  pp.  283-289. 

10.  Oyibo,  G.A.,  "The  Unified  Panel  Flutter  Theory  and  the  Damping 
Effects",  AFDC  paper  No.  82-002,  presented  at  the  Aerospace 


8 


11.  Weisshaar,  T.A.,  and  Ryan,  R.J.,  "Control  of  Aeroelastic  In¬ 
stabilities  Through  Elastic  Cross-Coupling”,  AIAA  paper. 

No.  84-0985,  May  1984. 

12.  Jensen,  D.W.,  and  Crawley,  E.F.,  "Frequency  Determination 
Techniques  for  Cantilevered  Plates  with  Bending-Torsion  Coup¬ 
ling",  AIAA  Journal,  Vol.  22,  No.  3,  March  1984. 


9 


4.0  RESEARCH  OBJECTIVES 


The  overall  goal  of  this  ohase  of  the  research  nroqram  is  twofold: 
first,  is  to  formulate  closed-form  solutions  to  the  aeroelastic 
divergence  and  flutter  of  aircraft  wings  with  warning  contraint 
and  elastic  coupling,  and  use  these  to  study  the  effects  of  con¬ 
strained  warnina  on  aeroelastic  response;  second  is  to  develop  new 
plate-beam  finite  elements  and  use  in  modal  analysis  of  vibration 
and  flutter  to  provide  an  independent  check  on  the  closed-form  so¬ 
lutions  results. 

The  results  of  the  efforts  towards  achieving  this  overall  objective 
are  the  subject  of  this  final  technical  renort.  Initially  ef¬ 
forts  were  concentrated  at  conducting  investigations  to  examine 
various  methods  (including  the  one  already  used  for  free  vibration 
studies)  with  which  to  formulate  the  closed-form  solutions  to  the 
divergence  problem.  Two  analytical  methods  were  thoroughly  exam¬ 
ined  to  determine  their  suitability  for  formulating  the  closed- 
form  solutions.  These  are:  (a)  the  method  of  "elimination", 

which  was  used  to  formulate  the  closed-form  solutions  to  the  free 
vibration  nroblem  durinq  the  first  phase  of  this  research,  and 
(b)  the  method  of  Laplace  Transformations.  The  latter  method  ap¬ 
pears  to  invoke  relatively  less  algebra,  and  so  was  implemented 
for  formulating  the  close-form  solutions  during  this  first  cruar- 
ter.  If  correct,  these  expressions  may  be  the  first  closed-form 
solutions  ever  derived  for  the  divergence  of  composite  aircraft 
wines  having  warping  constraint  and  elastic  coupling.  In  oarallel 
with  the  analytical  derivation,  the  numerical  approach  used  durinq 
the  first  phase  of  the  study  for  extracting  the  eigenvalues  from 
the  closed-form  expressions  for  the  free  vibration,  have  been 
carefully  studied  to  determine  how  it  (numerical  approach)  may  be 


used  for  extractinq  the  eigenvalues  (and  hence  flight  divergence 
sneeds)  from  the  closed- form  expressions  just  derived.  The  re¬ 
sults  extracted  were  found  to  be  in  agreement  with  previously 
published  results.  This  provided  the  necessary  confidence  in 
the  new  trends  predicted  by  the  new  theory.  Also  in  parallel 
is  the  development  of  a  new  plate  beam  finite  elements,  with 
which  to  independently  check  the  closed  form  solutions  being  gen 


erated . 


5 . 0  STATUS  OF  RESEARCH  EFFORTS 
5.1  INTRODUCTION 

During  the  year,  the  research  oroqram  progressed  accordina  to  olan 
resultincr  in  the  accomplishment  of  the.  goals  defined  for  the  year. 
Those  goals,  which  have  been  defined  in  Section  4  of  this  renort, 
basically  include  a  preliminary  formulation  and  investigation  of 
closed-form  solutions  to  the  aeroelastic  divergence  problems  of 
suoermanuverable  type  aircraft  wings  (fabricated  of  composite 
materials)  in  order  to  determine  an  effective  and  efficient  means 
with  which  to  accurately  estimate  and  study  the  effects  of  the 
"restraint  of  warping"  on  their  mechanism  of  aeroelastic  instabil¬ 
ity  . 

In  the  first  cruarter,  during  the  preliminary  study,  Prof.  G.  A. 
Oyibo  investigated  and  eventually  formulated  the  expressions  which 
seem  to  reoresent  the  closed  form  solutions  to  the  divergence  prob¬ 
lems  of  aircraft  wings  having  warping  constraint  and  elastic  coup¬ 
ling.  Prof.  J.  Bentson  explored  the  numerical  approach  used  for 
the  vibration  problem  in  order  to  determine  how  it  (numerical  ap¬ 
proach)  may  be  used  to  extract  the  divergence  eigenvalues  from  the 
closed  form  expressions  just  derived.  Prof.  T.  A.  Weisshaar  is 
also  continuing  his  investigations  into  how  a  new  beam-plate  finite 
elements  may  be  formulated  to  capture  the  constraint  of  warping  ef¬ 
fect  in  order  to  provide  results  which  may  be  used  to  independent¬ 
ly  check  the  closed-form  solutions  which  have  been,  and  are  being 
generated . 
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In  the  second  quarter  significant  programming  work  was  accom¬ 
plished  on  the  divergence  problem  and  this  effort  is  still  contin¬ 
uing.  Meanwhile  exciting  results  were  obtained  in  the  part  of  the 
research  work  that  calls  for  an  independent  verification  of  the 
analytical  closed  form  results  obtained  "thus  far  in  this  research 
program.  Therefore  it  is  a  pleasure  to  report  the  following: 

(i)  independent  results  obtained  from  MIT  seem  to  accurate¬ 
ly  verify  our  closed  form  results  as  is  shown  in  Fig¬ 
ure  6  in  the  copy  of  the  paper,  shown  below,  that  has 
been  sent  to  the  AIAA  Journal  to  be  considered  for 
publication. 

(ii)  an  elaborate,  very  careful,  independent  study  by  Prof. 

T.  A.  Weisshaar  and  his  team  at  Purdue  University 
(which  is  included  in  this  report  as  their  contribu¬ 
tion  to  the  program)  has  accurately  verified  our  closed 
form  results  and  some  of  their  implications.  For  ex¬ 
ample,  it  verified  our  proposed  physical  explanation 
(perhaps  the  very  first  explanation)  to  the  phenomenon 
of  "modal  transfer"  reported  previously  by  Purdue  and 
MIT. 

In  the  third  quarter,  more  exciting  conclusions  began  to- 
unfold,  that  might  potentially  change  the  way  dynamic is ts  look  at 
dynamics:  it  is  normally  thought  that  the  usual  way  one  gets  a 

complex  determinant  when  a  vibration  eigenvalue  problem  is  formu¬ 
lated  is  when  one  has  damping  in  the  system.  In  this  research  we 
have  shown  that  that  is  not  always  the  case.  We  have  now  seen 
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that  the  elastic  coupling  (D2g(or  D-^)  introduces  perhaps  some 
"conservation"  damping  characteristics  into  the  free  undamped 
vibration  of  an  anisotropic  wing.  This  results  from  the  com¬ 
plex  nature  of  the  determinant  and  the  consequential  fact  that 
the  bending  and  torsion  modes  are  out  of  phase  with  each  oth¬ 
er.  While  this  discovery  seems  to  puzzle  many  experienced  dy- 
namicists,  a  look  at  the  basic  equations  of  motion  would  in¬ 
dicate  that,  the  term  containing  (or  D,,)  is  an  odd  deriva- 

Zb  16  - 

tive  (similar  to  the  damping  terms).  Therefore,  it  would  seem 
that  this  physical  insight  which  is  perhaps  most  easily  extrac¬ 
ted  using  an  analytical  method,  is  not  only  meaningful  but  very 
important . 


In  the  final  ouarter,  the  divergence  problem  was  studied  more 
carefully  to  evolve  the  physical  insight  into  the  mechanism  of 
diveraence  instability  in  the  presence  of  warping  restraint  and 
elastic  coupling.  This  study  revealed  the  following: 


First,  it  is  seen  that  the  view  held  by  many  analysts  that  elast¬ 
ic  coupling  plays  a  significant  role  in  aeroelastic  tailoring  is 
verified . 

Second  it  is  seen  that  another  view  that,  higher  aeroelastic  sta¬ 
bility  boundaries  are  feasible  with  negative  elastic  coupling 
(than  positive  elastic  coupling),  is  basically  true,  but  up  to  a 
point.  It  is  further  seen  that  there  seems  to  be  a  limit  to  how 
negative  the  elastic  coupling  can  be  made  to  obtain  better  stabil¬ 
ity  boundaries  -  after  this  limit  a  further  negative  increment  in 
elastic  coupling  would  seem  to  result  in  lower  stability  bounda¬ 
ries. 
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Third  it  is  found  that  the  effective  aspect  ratio  defined  in  the 
first  phase  of  this  research  program  (>  )  for  simpler  models  can 
still  be  used  in  this  relatively  more  complex  model,  to  measure 
the  effect  of  warping  restraint  on  the  phenomenon  of  divergence 
instability.  The  results  show  that  ignoring  the  warping  restraint 
would  lead  to  conservative  estimates  for  the  divergence  instabil¬ 
ity  boundaries.  It  is  also  seen  that  the  restraint  of  warping 
effects  are  more  important  for  small  effective  aspect  ratios  and/ 
or  large  elastic  coupling. 
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5.2  ACCURATE  DIVERGENCE  THEORY  FOR  COMPOSITE  SUPERMANEUVERABLE 
AIRCRAFT  WINGS 
5.2.1  INTRODUCTION 

Modern  supermaneuverable  aircraft  concepts  benefit  a  great  deal 
from,  among  other  things,  significant  advances  in  materials 
technology  and  the  availability  of  more  accurate  aerodynamic  pre¬ 
diction  capabilities.  Supermaneuverability  as  a  design  goal  in¬ 
variably  calls  for  an  optimization  of  the  design  parameters. 
Optimization  may  be  partially  accomplished  for  example,  by  using 
composite  materials  to  minimize  weight.  Indeed,  it  has  been  known 
that  these  composite  materials  can  be  tailored  properly  to  resolve 
the  dynamic  or  static  instability  problems  of  these  types  of  air¬ 
craft.  The  concept  is  referred  to  as  aeroelastic  tailoring. 

While  aeroelastic  tailoring  has  tremendous  advantages  in  the  design 
of  an  aircraft,  the  analysis  which  provides  the  basis  for  the 
aeroelastic  tailoring  itself  is  generally  very  involved.  This  is 
rather  unfortunate  since  a  good  fundamental  physical  insight  of 
the  tailoring  mechanism  is  required  for  accurate  and  reliable  re¬ 
sults  . 

In  this  investigation  an  attempt  is  made  to  look  at  some  dynamics 
theories  that  can  be  used  to  understand  the  aeroelastic  tailoring 
mechanism.  Specifically,  the  accuracy  of  the  St.  Venant  torsion 
theory  which  is  relatively  simple  and  frequently  used  in  aeroelas¬ 
tic  analysis  is  examined  with  particular  reference  to  the  effects 
of  the  wings  aspect  ratio  as  well  as  other  design  parameters. 
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An  accurate  torsion/ twist  theory  is  particularly  significant  for 
supermaneuverable  aircraft  wings  since  supermaneuverability  is 
basically  characterized  by  high  angle  of  attack. 

Although  earlier  studies  (1,2,3)  have  indicated  that  the  St.  Ve- 
nant’s  torsion  theory  is  reasonably  accurate  except  for  aircraft 
wings  with  fairly  low  aspect  ratios,  the  theory  supporting  that 
conclusion  was  based  on  the  assumption  that  the  wing  is  construc¬ 
ted  of  isotropic  materials.  Basically  the  St.  Venant's  torsion 
theory  assumes  that  the  rate  of  change  of  the  wing's  twist  angle 
with  respect  to  the  spanwise  axis  is  constant.  This  assumption 
is  hardly  accurate  particularly  for  modern  aircraft  construction 
in  which  different  construction  materials  are  employed  and  the 
aerodynamic  loads  vary  significantly  along  the  wing's  span.  How¬ 
ever,  References  1-3  have  shown  that  (in  spite  of  such  an  inac¬ 
curate  assumption)  the  main  parameter  that  determines  the  accur¬ 
acy  of  the  St.  Venant's  theory  is  the  wing's  aspect  ratio.  Thus, 
it  was  determined  that  the  theory  is  fairly  accurate  for  moderate 
to  high  aspect  ratio  wings  constructed  of  isotropic  materials. 

In  recent  studies  (4,5,6)  however,  it  has  been  shown  that  for 
wings  constructed  of  orthotropic  composite  materials,  the  con¬ 
clusion  of  References  1-3  need  to  be  modified.  Rather  than  using 
the  geometric  aspect  ratio  of  the  wing  to  determine  the  accuracy 
of  St.  Venant's  twist  theory,  it  was  suggested  that  a  generic 
stiffness  ratio  as  well  as  an  effective  aspect  ratio  which  consid¬ 
ers  the  wing's  geometry  and  the  ratio  of  the  principal  directional 
stiffness  should  be  considered  in  establishing  the  accuracy  of  St. 
Venant's  theory. 
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The  present  investigation  is  related  to  the  studies  that  were  ini¬ 
tiated  in  References  5  and  6.  In  this  study  the  first  task  was  to 
examine  the  role  of  coupling  (both  mass  and  elastic  coupling)  on 
the  accuracy  of  St.  Venant's  theory  applied  to  static  problems. 

It  was  discovered  that  coupling  plays  a  very  significant  role  on 
the  accuracy  of  St.  Venant's  twist  theory.  The  second  task  was  to 
investigate  the  torsional  vibration  for  a  flat  plate  model  of  an 
aircraft  wing  fabricated  of  composite  materials  in  which  the  con¬ 
strained  warping  phenomenon  is  more  realistically  represented, 
with  particular  emphasis  on  higher  frequencies  and  to  compare  re¬ 
sults  with  those  from  a  representation  based  on  St.  Venant's 
theory . 

5.2.2  FORMULATION 

Consider  an  aircraft  wing  fabricated  of  composite  materials  and 
mathematically  idealized  as  a  cantilevered  plate  subjected  to  an 
aerodynamic  flow  over  its  surfaces.  The  mathematical  statement  of 
the  virtual  work  theorem  for  such  a  plate  model  is  well  known  and 
documented.  It  is  also  known  that  such  mathematical  statements 
of  the  virtual  work  theorem  for  a  laminated  plate  model  are  char¬ 
acterized  with  the  existence  of  so  many  variables  (in  the  state¬ 
ment)  ,  reflecting  the  various  directional  properties  for  the  lamin¬ 
ated  plate  model,  which  would  tend  to  interfere  with  any  physical 
insight  that  might  be  desired  from  a  phenomenological  analysis 
employing  such  a  mathematical  statement.  The  newly  discovered  af¬ 
fine  transformation  concept  (5,6  and  7)  was  developed  principally 
to  resolve  such  a  problem. 
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This  new  concept  therefore  can  be  used  to  evolve  the  mathematical 
statement  of  the  virtual  work  theorem  in  an  affine  space  given  by 
the  following  equation. 
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are  the  elastic  constants,  o,  is  the  material  density,  Ap 
is  the  differential  pressure  distribution,  w  is  the  displacement, 
t  is  the  time,  A  integrals  represent  area  integrals  and  h  is  the 
wing  box  depth. 

Equations  (1)  and  (2)  therefore  form  the  basis  of  the  newly  de¬ 
veloped  methodology.  The  equations  of  motion  of  a  plate  model  of 
an  aircraft  wing  can  now  be  derived  by  prescribing  a  realistic 
wing  displacement  and  using  Equation  (1). 

When  Equation  (1)  is  compared  to  its  physical  space  counterpart, 
it  is  seen  that  Equation  (1)  has  fewer  variables.  It  is  also  seen 
tha-!.  Equation  (1)  contains  only  non-dimensionalized  stiffness  quan¬ 
tities  (compared  to  dimensional  stiffness  quantities  in  its  physi¬ 
cal  space  counterpart) .  Another  feature  of  this  new  methodology 
which  makes  it  unique  is  that  the  non-dimensionalization  (a  conse¬ 
quence  of  the  affine  transformation)  is  accomplished  before  assum¬ 
ing  the  wing  deformations.  This  means  that  the  non-dimensionaliza¬ 
tion  is  independent  of  how  the  wing  deforms.  A  non-dimensionaliza¬ 
tion  scheme  that  depends  on  a  particular  assumption  of  the  wings 
deformations  could  lead  the  analyst  to  an  incorrect  physical  inter¬ 
pretation  of  results,  since  the  wing's  deformations  assumptions 
have  inherent  errors  because  they  are  based  on  the  analyst's  judg¬ 
ments  and  experience.  This  observation  may  become  clearer  during 
the  evolution  of  a  warping  parameter  with  which  to  study  the  ef¬ 
fects  of  the  warping  constraint  phenomenon  on  the  status  and  dyna¬ 
mics  of  a  wing  fabricated  of  composite  materials  later  in  this  study. 
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If  the  chordwise  curvature  is  neglected  in  an  initial  approxima¬ 
tion,  the  wing's  deformation  may  be  assumed  as  follows: 
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where  hQ  and  are  the  bending  and  twisting  displacements,  re¬ 
spectively  . 


It  can  be  shown  that  when  Equation  (3)  is  substituted  into  Equation 
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and  ,  the  following  equations  of  motion  are  obtained: 
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5.2.3  EVOLUTION  OF  WARPING  PARAMETERS 

The  evolution  of  the  warping  parameter  with  which  to  study  the  con¬ 
strained  aeroelastic  warping  phenomenon  for  wings  fabricated  of 
composite  materials  is  a  process  that  depends  on  the  sophistication 
of  the  wing's  mathematical  model;  whether  coupling  effects  are  in¬ 
cluded,  whether  the  wing's  chordwise  curvatures  are  included  and 
so  on.  Therefore,  any  warping  parameter  is  as  good  as  the  corres¬ 
ponding  wing's  displacements  assumptions.  However,  Equation  (1) 
makes  it  possible  for  the  analyst  to  determine  its  effective  inde¬ 
pendent  variables  even  before  the  displacement  assumptions  are  made. 

By  non-dimensionalizing  the  spanwise  space  variable  in  Equation  (4) , 
depending  on  whether  one  is  interested  in  the  static,  dynamic, 
coupled  or  uncoupled  displacements,  one  of  the  following  warping 
parameters  may  be  useful. 
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where : 


-  D  ( 1  - <  )  ;  ^2  =  2  (10) 

'  co 

(?0/cQ)  is  defined  as  the  wing's  effective  aspect  ratio  and  Dq* 
and  L  are  the  generalized  stiffness  and  coupling  ratios,  respec¬ 
tively  (defined  in  earlier  work  such  as  References  5  and  6) . 

Equations  (7)  thru  (9)  represent  the  appropriate  warping  parame¬ 
ter  for  dynamic  deformation,  static  displacement  with  elastic 
cross-coupling,  and  static  deformation  with  "geometric"  coupling 
(e  1  1) ,  respectively. 

It  was  discovered  in  this  study  that  evolving  the  warping  parame¬ 
ter  in  a  manner  shown  in  Equations  (1)  thru  (3),  should  enable 
one  to  investigate  the  effects  of  warping  on  the  composite  wing's 
dynamics  (or  the  accuracy  of  St.  Venant's  theory)  effectively. 
From  the  lamination  theory  for  composites  it  is  known  that  while 
Dq*  and  (£q/co)  are  always  positive,  L  and  can  be  positive  or 
negative.  However,  from  Equations  (1)  and  (3),  it  is  clear  that 
whether  a  composite  wing  has  positive  or  negative  coupling,  the 
warping  effect  (in  terms  of  X  )  is  unchanged. 
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Using  the  Laplace  transform  method  on  the  divergence  form  of  equa¬ 
tions  4  in  which  strip  aerodynamic  theory  is  used,  the  following 
Laplace  function  becomes  important  in  the  divergence  problem. 
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P0425  +  lP0426 


TQ4(~b2+ia2) 
F012{"b2+ia2)  '  042 
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P0425  +  (a2}  P0426 
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P043  P0435  +  (a3J  P0436 
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I 
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or  in  proper  notation 


or 


F2i(s) 


217 


P 


214 


215 


s+8 


14 


s+8 


15 


Z 


f21nS  +  P21n 
(s+Sn)2+a„2 


n=l ,2,3,6 


(28) 


fmn(s) 


MN7 


MN4 
s  +  8  , 


MN5 

s+8c 


Z 
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MNh  MNh 
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n=l ,2,3,6 


(29) 


MN  =  21,22,23,24 
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F21n(s) 


F22n(s) 


F23n(s) 


F.,  (s) 
24n 


(30) 


F217(S)  =  F210(S)/S  '  F214(S)  =  F210 (s) 7 (s+64) 

(31) 

F215^  =  F210^S^S  +  B5^ 


30 


(32) 


T2i(°) 


t21  (-e4) 


P217  "  '  P214  =  F214(-84) 
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'215  =  ^lS^V 
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or 
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21n5  x*21n6  F 
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n 
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(35) 
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p245  -  F215(-B5) 


24n5 


T24(-Bn+ian) 
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,  2  b4 

b4  P014e 


2 
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om 
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om 
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om 


or 
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3  b4 

F  (1)  -  bfp  .e 
ora  4^om4 


m  =  1 , 2 , 3 , 4 
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* MNn  MNn  n 
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2 


,  -3 

.3  m 
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n=l , 2 , 3 , 6 
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F11  F21  (1)  ;  F12  F22  ^  ;  F13  F23  (1)  '  F14  ~  F24(1^ 


F21  F01(1)  '  F22  F02(1)  ;  F23  F03(1)  ?  F24  _F04(1) 


alF21(1^  +  a5F01(1)  '  F32  alF22 “a5F02 ; 


F33  alF23 ^1^“a5F03 ^ 


F34  alF24(1^  +  a5F04(1) 


a  "2  _  a  -2 
(a4  _  )  F01(l)-a3F  (1)  ;  F 42  =  a3FQ2 ( 1) - ( a4~=  )F02(1) 

al  al 


a  2  a  2 

F43  =  a3F03(1)  ^ a4~—  *  F03^;  F44  =  ^a4~-  ^ F04 ^ ^ “a3F04 ( ^ 

al  al 


Determine  eigenvalues  that  makes  the  following  determinant  zero 


11 

F12 

F13 

F14 

21 

F22 

F2  3 

F2  4 

31 

F3  2 

F3  3 

F3  4 

41 

F42 

F4  3 

F  4  4 

Equation  50  is  therefore  expected  to  represent  the  closed-form  ex¬ 
pression  for  the  divergence  eigenvalues.  Numerical  methods  are  now 
being  used  to  e^-txact  these  eiaenvalues  in  order  to  obtain  the 
divergence  sdcoJs. 
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RESULTS 


The  results  obtained  so  far  seem  to  verify  existing  results  and 
establish  new  trends.  Due  to  the  fact  that  the  evolution  of 
these  new  trends  are  yet  to  be  fully  completed,  only  a  summary 
of  these  results  are  Dresented  below. 

The  results  obtained  by  solving  equation  50  were  carefully  to 
evolve  the  nhysical  insight  into  the  mechanism  of  divergence  in¬ 
stability  in  the  presence  of  warping  restraint  and  elastic  coun¬ 
ting.  This  study  revealed  the  following,  as  can  be  seen  from 
fiaures  1-6. 

First,  it  is  seen  that  the  view  held  by  many  analysts  that  elas¬ 
tic  couplinq  olays  a  significant  role  in  aeroelastic  divergence 
tailorino  is  verified. 

Second,  it  is  seen  that  another  view  that,  higher  aeroelastic 
divergence  stability  boundaries  are  feasible  with  negative  elas¬ 
tic  couplina  (than  positive  elastic  coupling) ,  is  basically  true, 
but  uo  to  a  ooint.  It  is  further  seen  that  there  seems  to  be  a 
limit  to  how  neoative  the  elastic  coupling  can  be  made  to  obtain 
better  stability  boundaries  -  after  such  a  limit,  a  further  nega¬ 
tive  increment  of  elastic  couoling  would  seem  to  lower  the  stabil- 
i ty  bounda  r i e  s . 

Third,  it  is  found  that  the  effective  aspect  ratio  defined  in  the 
first  ohase  of  this  research  program  {)  ) ,  for  simoler  models  can 
still  bo  used  in  this  relatively  more  complex  model,  to  measure 
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the  effect  of  warping  restraint  on  the  phenomenon  of  divergence 
instability.  The  results  show  that  ignoring  the  warpinq  re¬ 
straint  would  lead  to  conservative  estimates  for  the  divergence 
instability  boundaries.  It  is  also  seen  that  the  restraint  of 
warping  effects  are  more  significant  for  small  effective  aspect 
ratio  ('  )  and/or  large  elastic  coupling. 
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L 

Figure  6  Variation  of  divergence  dynamic  pressure  with  coupling 
for  large  effective  wing  aspect  ratio 
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5.3 

Exact  Solutions  to  Aeroelastic  Oscillations  of  Composite  Aircraft 

•  •  t 

Wings  with  Warping  Constraint  and  Elastic  Coupling 

*  * 

Gabriel  A.  Oyibo  and  James  Bentson 

Polytechnic  University 
Farmingdale,  New  York 

Abstract 

Exact  solutions  within  the  framework  of  standard  aeroelastic 
bending  and  twisting  assumptions  are  found  to  the  free  oscilla¬ 
tions  of  composite  aircraft  wings  having  warping  constraint  and 
elastic  coupling.  The  problem  is  treated  as  a  regular  boundary 
value  problem  consisting  of  two  fourth  order  partial  differen¬ 
tial  equations  coupled  by  the  presence  of  elastic  coupling. 

This  system,  which  is  linear,  therefore  is  equivalent  to  an 
eighthorder  ordinary  differential  equation  .  Classical  linear 
"operator"  method  is  therefore  used  to  extract  fundamental  solu¬ 
tions  which  are  superimposed  appropriately  to  obtain  an  exact 
functional  form  for  the  mode  shapes.  These  mode  shapes  are 
therefore  made  to  satisfy  the  necessary  boundary  conditions,  a 
process  that  leads  to  the  formulation  of  the  required  eigenvalue 
problem.  The  eigenvalues  are  extracted  numerically  by  using  ap¬ 
propriate  ordering  of  the  eight  roots  of  the  operator  equation. 

The  bending-torsion  frequencies  obtained  as  a  result  of  this  an¬ 
alysis  are  compared  favorably  with  existing  results.  New  in¬ 
sights  made  possible  by  these  results  which  are  preliminary,  ap- 

~ Research  sponsored  by  the  Air  Force  Office  of  Scientific  Research 
(AFOSR) ,  under  Contracts  F49620-85-C-0090  and  F49620-87-C-0046 . 

Associate  Professor  Dept,  of  Aerospace  Engineering 
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pear  to  be  that  (a;  the  first  coupled  frequency  decreases  with 
increasing  coupling  and  (b)  the  phenomenon  of  modal  transforma¬ 
tions  found  by  earlier  investigators  is  explainable  in  terms  of 
some  conservative  inter-modal  energy  transfer. 
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5.3.1  NOMENCLATURE 


c  ,  c 
o'  o 


=  chordwise  integrals 

=  affine  space  half-chord  and  chord,  respectively 
=  elastic  constants 

=  parameter  that  measures  the  location  of  the 
reference  axis  relative  to  mid-chord 


El,  GJ 


=  bending  and  torsional  stiffness,  respectively 


=  wing  box  depth 


(ho,cto) 


=  affine  space  bending  and  torsional  displacement, 
respectively 


K ,  So 


=  elastic  coupling  and  warping  stiffness,  respec¬ 
tively 


=  elemental  stiffness  parameter 


=  Strouhal  number 


L  .M„ 
o  o 


affine  space  running  aerodynamic  lift  and  mo¬ 
ments,  respectively 


=  affine  space  half-span  for  the  wing 
=  affine  space  mass  per  unit  span 


(Ap, Apo) 


=  differential  aerodynamic  pressure  distributions 
in  physical  and  affine  space,  respectively 


=  time 

=  virtual  work  expressions  in  physical  and  affine 
space,  respectively 

=  flutter  speed 


(x, y , 2 , ) ,  (xo,yo,2Q) 


=  physical  and  affine  space  coordinates,  respec¬ 
tively 


r , L.  , , D* , D*  , 
12  O 


=  displacement  shape  functions 

=  generic  nondimensionalized  stiffness  parameters 
=  affine  mass  ratio  parameter 


uij,e 


=  Poisson  ratios  and  generalized  Poisson's  ratio, 
respectively 


=  divergence  parameter 


5.3.2 

Introduc  tion 


Perhaps  one  of  the  more  elusive  aspects  of  supermaneuverabil¬ 
ity  as  a  design  concept  is  its  aeroelastic  implications.  One  gen¬ 
erally  accepted  definition  of  supermaneuverable  aircraft  is  that  it 
is  designed  to  operate  at  high  angles  of  attack.  Strictly  speaking 
high  angle  of  attack  problems  are  nonlinear.  However  due  to  the 
high  degree  of  complexities  involved  in  dealing  with  nonlinear  aero¬ 
elastic  problems, an  average  aeroelastician  would  prefer  to  deal  with 
a  linearized  version  of  the  problem  (at  least  as  a  first  approxima¬ 
tion)  .  If  linear  aeroelastic  equations  are  used  under  such  condi¬ 
tions,  at  least  it  should  be  assumed  that  the  high  angle  of  attack 
would  introduce  large  twisting  displacements  which  would  imply  that 
terms  containing  twisting  displacement  should  be  retained.  Even 

under  low  angle  of  attack  assumptions,  the  early  works  of  Reissner 

1  2 
and  Stein  and  later  works  of  Libescu  et  al  have  shown  that  for 

metal  wings  there  are  conditions  under  which  the  so  called  St.  Ve- 
nant's  torsion  principle  is  inapplicable.  This  is  when  the  res¬ 
traint  of  warping  effect  is  important  and  a  more  accurate  analysis 
would  need  to  include  a  higher  order  term  involving  the  twisting 
displacement.  Although  the  retention  of  such  a  term  implies  solv¬ 
ing  fourth  order  (instead  of  second  order)  differential  equations 
for  the  twisting  and  bending  displacements,  the  equations  can  be 
easily  decoupled  for  metal  wings.  However  for  composite  wings,  the 

decoupling  of  these  equations  is  neither  easy  practically  nor  is  it 

3-7 

even  desirable  from  aeroelastic  tailoring  standpoint.  These  stu¬ 

dies  have  also  shown  that  the  restraint  of  warping  is  very  important 
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in  composite  winas .  Tnerefore  i  t  would  seem  that  an  aeroelastic 
analysis  of  a  supermaneuverable  (hinh  angle  of  attack)  aircraft 
wing  fabricated  of  composite  materials  would  need  to  consider 
the  effects  of  restraint  of  warping  as  well  as  elastic  coup¬ 
ling  . 

Previous  investigation  of  this  latter  problem  (free  vibra- 
3  4 

tion)  at  MIT  '  used  analytical  methods  to  solve  the  decoupled 
problem,  while  numerical  methods  were  utilized  to  solve  the  cou¬ 
pled  problems.  Consequently  general  results  were  presented  for 
the  decoupled  problem  while  representative  results  were  presen¬ 
ted  for  the  coupled  problem. 

In  this  paper  the  coupled  free  vibration  is  treated  analyti¬ 
cally  as  a  pair  of  coupled  fourth  order  (differential  equations) 
boundary  value  problem  to  which  exact  closed  form  eigen-solutions 
are  sought.  The  enforcement  of  the  necessary  boundary  conditions 
resulted  in  a  fairly  complicated  transcendental  function  to  be 
used  to  determine  the  required  eingenvalues  from  which  the  natur¬ 
al  frequencies  are  to  be  obtained.  This  transcendental  function 
was  complex  in  contrast  to  its  decoupled  counterpart  (which  is 
real) .  That  should  be  indicating  the  presence  of  the  phase  angle 
that  exists  between  the  twisting  and  bending  displacements.  A  com¬ 
parison  with  a  damped  (decoupled)  system  in  which  complex  determin¬ 
ant  signifies  phase  angles  between  damping  and  other  forces,  led  us 
to  the  formulation  of  an  explanation  for  the  "modal  transformation" 
phenomenon  which  was  reported  in  studies  at  MITJ'4  and  Purdue  (which 
seemed  to  have  lacked  explanation  until  this  study) .  The  explana¬ 
tion  is  that  the  modal  transformation  may  be  viewed  as  a 
form  of  steady  state  conservative  (energy  stays  in  this  system 
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since  there  is  no  clampinq)  inter- 


modal  energy  transfer  between  the  vibration  modes.  In  fact 

g 

work  currently  in  progress  at  Purdue  seem  to  support  and  con¬ 
firm  this  explanation.  The  results  which  favorably  compared 

3  4 

with  those  obtained  at  MIT,  '  also  revealed  that  coupling  has 
a  tendency  to  lower  the  first  coupled  natural  frequency  of  a 
composite  aircraft  wing.  In  fact  it  is  seen  that  a  substantial 
amount  of  coupling  could  reduce  the  first  coupled  frequency  to 
almost  zero  (hence  a  possibility  of  coupling  with  rigid  body 

modes)  . 

5.3.3 

Problem  Statement 

For  a  composite  aircraft  wing  cantilevered  at  the  root  as 
shown  in  Figure  1,  the  virtual  work  theorem  in  the  physical 
space  is  given  by 


6U  =  0  =  j 


+  d22( 


fill 

0  A  L 


Dll(w'xx)  +  2D12  W'xx  W,yy 


+  4D16  w'xx  w'xy  +  4D26  w'yy 


w,  )  2 

yy 

4D66  (w'x  )2  dxd^dt  ~  if  f  I  p  h  W2dxdydt 

t  •'o  •'a*' 

•/// 


(1) 


+  ADCC.  (w, 

Ap6w  dxdydt 


where : 


D  are  the  elastic  constants,  p,  is  the  material  density,  Ap 
is  the  differential  pressure  distribution,  w  is  the  displacement, 
t  is  the  time,  A  integrals  represent  area  integrals  and  h  is  the 
wing  box  depth. 

Using  the  following  affine  transformation  of  variables; 


D1 1  !/4 

x  =  (-=i)  x 
D22 


y  =  y. 


z  =  z  (2a) 
o 


s  1 


then  in  affine  space  the  virtual  work  theorem  becomes: 


where 


6U  =0=2 
o 


III 

0Y 


(w,x  )-  +  2D* 

o  o 


(1"E) (W'x  y  ) 
oJ  o 


+  eW'v  v  W '  \r  W  +  (W*„  )  *"  +  L,W,  W  + 

o  o  yoyo  yoyo  ^  xoxo  Xoyo 


+  L2W'y  y  %y  [  dxodyodt 
J  oJ  o  oJ  o 


(2b) 


II! 


iff 


dxodyodt  +l  II  Ap^6wdxody^dt 


O  O  J  O 


1/4 


u  =  JL.  (^22 

°  D22  \D11 


D 


;  eD 


*  _ 


12 


(DUD22) 


4D 


L„  = 


26 


D*  = 


P12  +  2D66 

(D11D22)1/2 
4D 


1/2 


L1  = 


(D11)3/4(D22)1/4 


(D  )1/4(D  )3/4 

^Uir  'U22 


= 


D22 


P  = 

°  D22 


(3) 


If  the  affine  space  equivalent  of  the  standard  aeroelastic 
displacement  assumptions  is  made,  i.e.. 


w(xo'yo'fc)  = 


ho{yo't}  +  xoao(yo't) 


(4) 
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where  hQ  and  .1  are  the  bendinq  and  twistinq  displacements  re¬ 
spectively,  then  it  can  be  shown  through  the  use  of  the  calcu¬ 
lus  of  variations  that  a  coupled  set  of  aeroelastic  equations 
of  motion  for  the  composite  aircraft  wing,  in  which  the  re¬ 
straint  of  warping  and  elastic  coupling  effects  are  accounted 
for  is  given  by 


a.h  1v  +  a  a  iv  -  a.  a’n  +  p  a.h  +  p  a.ci  «  L 
1  o  2  o  >>  0  o  1  o  2  o  o 


a2hoV+65holi+  +  P°33“o  +  Po a2’h’ 


(5) 


_  -  M 

o  o 


at  y  =0 
o 


with  boundary  conditions 


s 

o 

II 

r: 

h'  = 

II 

1 

s 

o 

o 

.  o 

o 

at  y  = 
J  o 

a,«  + 

£ 

o 

a0h"  - 

■  L0a0a^ 

3  o 

2  o 

2  2  o 

a„a  + 

a,h~ 

-  aca 

2  o 

1  o 

5  o 

a«h  + 

a  ->ct 

"  +  a,h' " 

2  o 

3  o 

5  o 

l-i  U  Cl  r-  Cl 

1  o  5  o 


2  o 


l  ,  u. 

4  o 


(6) 
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where : 


a  , 


f  dx 

Jec.  0 


■Jj' 


o 

o  o 


J  tc 


X  2dx 
0  0 


3  4  *  2 


f  f°  D* 


(l-€)dxf 


.  p 

^ec  c 


.  /p  dx  .  a, 

o  ^ec0  o  o  i  5 


•  f  c°  4dx< 

J  ec 


(7) 


f: 

ec. 


x0  A  podx0 


-  ®  <  e  <  o 


c0  =  1-e 


1  1  sy0-  1  1  *  aT 


For  free  vibrations , if  a2  (through  the  geometric  construction 
of  the  wing)  is  made  to  be  zero,  equations  5  reduce  to 


IV 


a,h  -  aca  +  pa, li  =  L 
lo  5o  olo  o 


(8a) 


iv  ,  ,  ^  - 

a^a  +  achi 
3  o  5  o 


-a,  a'"  +pa0a  =  M 
4  3  o  o 


with  boundary  conditions 
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at  y 


0, 


ho  =  O'  ho  =  0'  aQ  =  0,  a'  =  0 


at 


l 

o 


a 


o 


>  >  y 

a,-a 
5  o 


=  0, 


h-*  >» 

X  o 


aj-a'  =  0 
o 


a a"  '  +  a 
3  c 


5ho 


a4ao 


0 


(8b) 


It  may  be  stated  here  that  the  restraint  of  warping  effect  is  rep¬ 

ly 

resented  by  the  product  of  a..  and  a3  while  the  elastic  coupling 

effect  is  represented  by  the  a^  terms. 

5.3.4 

Methods  of  Solution 

Two  methods  for  solving  equations  8  are  examined  in  this  pa¬ 
per.  These  are  (a)  an  exact  closed  form”  approach  and  (b)  a  "semi- 
exact  closed  form  approach.  The  exact  closed  form  approach  is  de¬ 
fined  here  as  one  in  which  explicit  expressions  are  derived  for 
the  eight  roots  of  the  eighth  order  operator  equations  represent¬ 
ing  equations  8a,  and  through  the  superiraposi tion  of  fundamental 
solutions  corresponding  to  each  of  the  eight  operator  roots,  the 
boundary  conditions  8b  are  satisfied.  The  semi-exact  closed  form  ap¬ 
proach  is  the  same  as  (a)  except  that  the  roots  of  the  operator 
equation  are  determined  numerically  through  the  usage  of  some 
standard  root  extraction  subroutines. 

In  either  case,  to  solve  for  the  operator  roots  of  equation 
8a,  it  may  be  rewritten  in  operator  form  as  follows, 
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0 


(9a) 


I*  -i  h  Xl  r-|  IX  — 
1  o  2  o 


L^h  +  L0a  =  0 
2  o  3  o 


where  the  operators  Ln  (n  =  1,2,3)  are  given  by 


4  3 

7  _  9  2  7 

L1  “  al  4  ”  alw  po  L2  "  a5  — 3 


L  -  a  34  -  a  ^ 

j  ~  a3  4  a4  2 


a,u5  p 

3  Ko 


and  h  =  h 
o  o 


—  iwt 
a  =  a  e 
o  o 


ai  is  vibration  frequency 


Following  the  method  used  in  Hilderbrand,  if  and  L2  are  com¬ 
mutative  (i.e.,  L^L2  =  L2Ll^  or  ^2  an<^  L3  are  conunutative  (i.e., 
L2l3  =  L3L2)  either  hQ  or  aQ  may  be  eliminated  and  equations  9  is 
reduced  to  one  eighth  order  operator  equation  in  either  aQ  or  hQ, 
respectively.  This  can  be  accomplished  by  using  L2  to  operate  on 
9(a)  and  to  operate  on  9(b)  and  subtracting  9(a)  from  9(b)  or 

by  using  L3  to  operate  on  9(a)  and  L2  to  operate  on  9(b)  and  sub¬ 
tracting  9(a)  from  9(b).  The  result  of  such  an  operation  in  which 
hQ  is  eliminated  is  given  by 
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0 


(x  +  F,x  +  F_x  +  F.x  +  F  )  a  = 

j  2.  loo 

2 

where,  F^  =  -16Ac^,  F2  =  -k  ,  F^  =  8  (kA  )  ^ 


F  =  7—  ,  X  =  (7 - ) 

o  4  3y 

1  o 


It  should  be  noted  here  that  because  the  method  of  elimination  of 
hQ  is  through  a  differential  operator,  superfluous  solutions 
are  to  be  expected  (e.g.,  operator  equation  11  is  eighth  order 
instead  of  fourth  order) .  These  superfluous  solutions  are 
eliminated  by  enforcing  some  consistency  conditions  on  the  so¬ 
lutions  as  described  below.  Based  on  equation  11,  an  appropriate 
form  for  the  exact-closed  form  (or  complete)  solution  for  the 
mode  shapes  of  the  free  vibration  problem  proposed  in  this  paper, 
which  can  also  be  made  to  satisfy  the  necessary  boundary  condi¬ 
tions  is  given  by. 


hQ  =  coshB1yQ  +  A2  sinhB^yQ  +  A3  cos82yo  + 
A4  sin82Y0  +  A5  coshB3yQ  +  A&  sinh83y0+ 
A?  cos84y0  +  Ag  sinB4yQ 

aQ  =  cosB1y0  +  B2  sinhB1yc  +  B3  cosB2y0  + 
B4  sinB2yQ  +  B5  coshB3yQ  +  Bg  sinhB3yQ  + 

B7  cos84yo  +  B8  sinB4yo 


(12(a)) 


(12(b)) 


where 
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B  =  x  1/2 
n  n 


n  =  1,2, 3, 4 


and  xn  can  be  obtained  either  from  the  numerical  solution  of 
equation  11  or  according  to  the  method  described  in  Abramo- 
witz  and  Stegun,10  as  follows 
Define 


X_  =  U  /c  )  /Id 

c  o  o  2  o 


Xc  -  (VCo’  v4(Do  “  2  )  (14) 


_  2 

4  (X^  +  X  )— 8/3X  X2  +  x-y 
c  c  c  c  2  / 


1  2  lc 

i  (32J/V  +  £  ) 

O  C  C  j 


Sx  =  2k  kP-L  +  (k2P12  - 


S2  =  2k 


8P23) 


fkpr 


(k2Px2  -  8P3) 


1/2 | 1/3 


ui  —  sx  +  s2  -  2 


u2  “  "  2  (S1+S2)_  ^  +  1  2  (S1-S2) 


u3  “  "  I  {S!+S2)  “3  1  2  (S1-S2) 
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f 


3,4 


f 


5,6 


+  u  . 
1 


2 


(18) 


x 


1,2 


x 

3,4 


(19) 


where  (i  =  1,2,3) 


is  the  root  that  makes  f,  .  c  .  all  real. 

3,4 ,5,6 


5.3.5 

Consistency  Conditions 

From  equations  12(a)  and  12(b)  it  is  seen  that  there  are  six¬ 
teen  arbitrary  integration  constants  as  opposed  to  the  expected 
eight  constants.  The  additional  eight  constants  have  been  intro¬ 
duced  superfluously  as  a  result  of  the  differential  operation 
which  was  done  in  order  to  eliminate  one  of  the  dependent  varia¬ 
bles  in  the  two  coupled  differential  equations  (equations  9) .  In 
order  to  get  rid  of  these  superfluous  solutions  it  is  necessary  to 
enforce  some  consistency  conditions.  This  may  be  accomplished  by 
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substituting  equations  12(a)  and  12(b)  into  one  of  the  equations 
9  and  requiring  the  equation  to  be  satisfied  identically.  This 
procedure  shall  establish  a  set  of  explicit  relationships  be¬ 


tween  the  constants  A  and  B  in  which  A  can  be  determined  in 

n  n  n 

terms  of  B  or  vice-versa, 
n 

When  the  superfluous  solutions  are  eliminated,  equations  12 
can  now  be  used  to  satisfy  the  boundary  conditions  for  this  prob¬ 
lem  (equations  8b) .  Consequently  the  condition  for  nontrivial  so¬ 
lutions  is  enforced  to  obtain  the  transcendental  functional  ex¬ 
pression  for  determining  the  eigenvalues  of  this  problem. 

5.3.6 

Eigenvalues 

The  following  steps  and  definitions  are  carried  out  in  order 
to  obtain  the  transcendental  functional  expression  for  determining 
the  necessary  eigenvalue  (for  this  eighth  order  boundary  value 
problem)  from  which  the  coupled  natural  frequencies  may  be  ob¬ 
tained 


Define 


(20) 


t 

n 


16Uc2-ic2i 


82  +  (-l)n+1  — 
n 


B 


I  [8n2  +  (-l)n16Xc2) 


n 


(21) 


fin 


F51  =  -613a15sinh81+623a35sin62+&33sinhB3 

3  3 

F61  =  B13a26cosh81-B2  a46cos82+33  cosh83 

7  3  3  • 

F?1  =  B1  a17sinhB1-62  a37sinS2+&4  sinB4 

F81  =  613a28coshB1-B23a4gcosB2-B43cosB4 
F52  =  h£al5sinhB1-h2a35sinB2-h3sinhB3 
F62  =  -h3a26coshB1+h2a46cosB2-h3coshB3 
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F 


72 


F 


82 


F 


53 


F 


63  ~ 


F 


73 


F 


83  ~ 


F 


54 


F 


64  " 


F 


74 


F 


84 


“hla17sinhBl+h2a37sini32-h4sin,34 
-illa28COSh^l"t’h2a4  8GOS^2+h4COS^4 


t ^ ^ ^ ^ co s 8 ^ t cosh 8 ^ 


tia26S;''nh^l-t2a46S^nl52+t3s;'"nll^3 


tla17COS^1^l+t2a3  7cos^2_  t4COS^4 


tla28sin^^ l_^2a4  8sin^2~^4sin^4 


tial5COS^^l+t2a3  5COS^2~t3COsIl^3 
-tla26Sanh®l+t2a46San^2-t3S:’Lnh^3 


- t^a^^cos8^-t2a2yCOs32+t^cosB^ 
~^la2  8s;''nll^l+:^2a4  6'3an^2  +  t4S'''n^4 


(24) 
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(25) 


a 
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F61F72  F 


F51F62-F 


6  2F7  1 
52F6 1 


a58 


F61F82  F62F81 
F51F62_F52F61 


a 
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F52F71  F 


F51F62-F 


51F7  2 
52F61 


68 


F52F81  F51F82 
F51F62-F52F61 


Then,  kn  are  given  by  the  roots  of 


F  = 


(F53a57+F63a67+F73) (F54a58+F64a68+F84) " 


(26) 


(F54a57+F64a67+F74}  (F53a58+F63a68+F83)  ° 


Equation  26  is  therefore  the  exact  closed  form  transcendental 
functional  expression  from  which  the  eigenvalues  kn  may  be  ex¬ 
tracted.  The  coupled  natural  frequencies  are  related  to  these 
eigenvalues  by  the  following  expression. 


U) 


n 


w 


(27) 


5.3.7 

Computations 

The  extraction  of  the  eigenvalues  from  the  exact-closed  form 
transcendental  expression  in  equation  26  proved  to  be  a  very  chal¬ 
lenging  computational  exercise  due  mainly  to  its  complex  nature, 
the  existence  of  branch  points,  the  necessity  to  order  the  opera¬ 
tor  roots  appropriately  and  the  existence  of  numerical  noise. 
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The  characteris tic  roots  for  the  equations  of  motion,  x  , 

n 

were  extracted  in  two  different  ways  in  order  to  assure  accuracy: 
one  method  was  by  using  a  Jenkins-Traub  computational  method  from 
the  IMSL  library.  The  second  method  was  by  using  an  exact-closed 
form  approach  outlined  in  Abramowitz  and  Stegun.^  When  these 
roots  are  being  extracted  numerically,  they  do  not  necessarily 
come  out  in  a  continuous  manner.  Therefore  a  subroutine  that  re¬ 
orders  them  so  that  they  become  continuous  with  respect  to  the 
eigenvalue  is  also  employed.  Once  this  reordering  of  these 
roots  is  completed,  the  parameters  needed  for  computing  the  trans¬ 
cendental  function,  F,  in  equation  26  are  computed.  Finally,  the 
transcendental  expression  itself  is  computed  and  the  roots,  k^, 
are  found  numerically. 

One  interesting  result  is  that  the  values  of  F  are  complex 
here.  It  is  known,  however,  that  each  of  the  two  uncoupled  prob¬ 
lems  (bending  or  torsion)  when  treated  separately  has  a  real 
transcendental  expression  for  extracting  the  eigenvalues.  This  be¬ 
havior  of  the  transcendental  expression  for  extracting  the  coupled 
eigenvalues  (i.e.,  being  complex  as  opposed  to  being  real)  was  the 
first  hint  that  led  the  first  author  to  examine  any  possible  math¬ 
ematical  similarity  between  the  problem  at  hand  (a  coupled  non- 
damped  oscillations  problem)  and  simple  damped  oscillations  prob¬ 
lem  in  which  the  expression  for  determining  the  eigenvalues  is  in 
general  complex  due  to  the  need  to  determine  the  oscillation  fre¬ 
quency  and  the  amount  of  damping  in  the  mode,  etc.  In  the  coupled 
problem  the  complex  nature  of  this  transcendental  expression  for 
extracting  the  eigenvalues  seems  to  be  basically  a  reflection  of 


64 


the  phase  that  normally  should  exist;  between  the  bending  and  tor¬ 
sional  modes.  Realizing  that  in  the  damped  case,  there  is  a  non¬ 
conservative  energy  transfer  between  the  oscillating  system  and 
its  environment  as  well  as  intermodal  energy  transfer,  it  became 
clear  to  the  first  author  that  in  the  undamped  coupled  system 
there  may  be  a  conservative  energy  transfer  in  the  oscillating 
system.  At  the  time  of  this  thought  it  seemed  to  the  first  auth¬ 
or  that  if  it  made  sense,  the  idea  should  provide  an  explanation 

for  the  modal  changes  which  were  noticed  during  the  studies  at 
3  4 

Purdue  and  MIT  '  as  coupling  in  the  system  was  changed.  This 

prediction  which  turns  out  to  be  useful  and  confirmed  by  other 

0 

recent  studies  from  Purdue  University  shall  be  discussed  in  more 
detail  below. 

The  most  convenient  way  to  obtain  the  eigenvalues,  kn ,  was 
found  to  be  by  graphical  means.  Thus  the  values  of  the  real  and 
imaginary  parts  of  the  complex  transcendental  function  F  are 
plotted  against  the  values  of  on  the  same  curve  as  shown  in 
Figure  2.  The  values  of  k^  at  which  the  real  and  imaginary  parts 
of  the  transcendental  function  are  simultaneously  equal  to  zero 


corresponds  to  the  desired  eigenvalues  for  this  coupled  problem. 

One  of  the  problems  with  the  computational  model  described 
above  is  that,  when  the  coupling  parameter,  ,  becomes  identical¬ 
ly  zero,  the  coupled  system  of  equations  becomes  computationally 
ill  conditioned  and  unsolvable.  To  circumvent  this  the  results 
for  the  uncoupled  case  Uc/Ac  =  D  are  obtained  from  a  series  of 
calculations  using  successively  srru  ’  ler  values  for  the  coupling 
parameter  L  .  This  led  us  to  accept  the  values  of  eigenvalues 
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for  zero  coupling  as  the  value  corresponding  to  the  limit  as  the 
coupling  approaches  zero.  Although  this  continuity  assumption 
seems  to  make  sense,  it  is  not  backed  by  a  rigorous  mathematical 
proof.  Luckily  it  was  possible  to  check  these  results  with  re¬ 
sults  generated  for  an  isotropic/metal  aluminum/wing  by  MIT,3,4 
and  the  agreement  was  found  to  be  good  for  the  cases  checked. 

5.3.8 

Results  and  Discussions 

The  natural  frequencies  w  for  the  coupled  bending- torsion 

oscillations  for  a  composite  aircraft  wing  in  the  presence  of 

elastic  coupling  and  warping  restraint  is  found  (as  shown  by 

1/2 

equation  27)  to  be  a  function  of  the  ratio  (D22/P)  <  the  length 

and  the  nondimensionalized  frequency  parameter  kR  for  the 
wing.  In  this  problem,  kn,  is  a  function  of  only  two  parameters, 
i.e.,  Ac  which  may  be  considered  as  an  effective  nondimensional¬ 
ized  aspect  ratio  and  Ac/Ac,  which  in  a  way,  measures  the  amount 
of  elastic  coupling  in  the  wing  (A  /A  =  1  for  zero  coupling) . 

It  is  therefore  seen  from  equation  27  that  in  order  to  increase 

10  one  needs  to  make  l  as  small  as  possible  and/or  make  (D_-/p) 
n  o  ‘‘ 

and  kn  as  large  as  possible.  Such  an  exercise  may  be  necessary 
when  a  tailoring  of  the  frequency  is  needed  to  avoid  instabili¬ 
ties  (e.g.,  very  low  structural  frequencies  may  provide  an  atmos¬ 
phere  for  a  coupling  between  the  flexible  modes  and  rigid  body 
motions,  which  in  turn  has  a  potential  to  result  in  instability). 
Tnis  kind  of  tailoring  is  made  convenient  through  the  use  of  equa¬ 
tion  27  in  which,  for  a  particular  wing  configuration  and  compos¬ 
ite  material,  every  variable  in  the  equation  shall  be  known  ex- 
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as  we  said  earli- 


cept  k  and  Obviously  if  we  want  hiah  w  , 

n  z  z  ^  n 

er '  ^22  be  made  to  be  as  high  as  possible  (and  of  course  p 

as  low  as  possible) .  Once  this  is  done  the  only  other  parameter 

to  be  tailored  is  k  . 

n 

The  plot  of  k^  as  a  function  of  Ac  and  Ac  is  shown  in  Figure 

3  for  all  configurations  having  low  camber  to  twisting  coupling, 

and  all  values  of  bending  to  twisting  coupling.  The  results  com- 
3  4 

puted  at  MIT,  '  which  were  used  to  verify  the  present  results, 

are  shown  in  Figure  3  as  well.  The  important  trend  made  visible 

is  this  investigation  as  shown  in  Figure  3,  is  that  k^  (and  hence 

u)jJ  decreases  with  increasing  nondimensionalized  coupling  , 

which  perhaps  may  be  a  more  effective  way  to  actually  measure  and 

compare  elastic  couplings,  D16  and  D2g .  For  example,  the  results 
3  4 

from  MIT  '  which  were  computed  for  some  representative  configura¬ 
tions,  in  dimensionalized  form  seem  to  represent  systems  with  a 
fairly  significant  variation  in  coupling  or  (depending  on 

the  coordinate  system).  However  when  nondimensionalized,  the  re¬ 
sults  as  shown  in  Figure  3  seem  to  show  little  variation  in  coup¬ 
ling.  In  fact  they  appear  to  be  so  close  to  the  zero  coupling 
case  (^c/^c  =  1)  or  isotropic  (or  metal)  case,  that  k^  for  a  metal 
or  isotropic  wing  should  be  a  good  approximation  (if  it  was  nec¬ 
essary  to  make  an  approximation) .  The  low  value  for  the  effective 
coupling  was  also  evident  in  the  nondimensionalized  results  from 
MIT3,4  where  the  bending  frequency  hardly  varied  with  material 
changes.  The  question  that  could  be  asked  is  therefore  "Do  all 
possible  composite  wing  configurations  result  in  very  low  effective 
nordimensionalized  coupling,  (A  /A  ~1) ?"  If  the  answer  is  "yes”, 

C  V/ 
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then  it  may  be  proposed  that  for  X  >.5,  k,  (A  , A  /A  )  may  be  ap- 
proximated  as  k^(Ac,l),  which  is  also  the  isotropic  or  metal 
value.  For  this  case  it  is  seen  that  the  computation  of  the  na¬ 
tural  frequencies  of  composite  aircraft  wing  having  aeroelastic 
oscillations,  merely  requires  the  computation  of  (D22/P)  for  a 

given  wing  half  span  l  ,  since  k  (A  , A  /A  )  is  approximately 

O  tl  o  c  c 

equal  to  k^  (Ac,l)  which  is  approximately  equal  to  a  constant 
(3.5)  for  n  =  1.  This  result  should  make  frequency  computations 
for  the  bending  mode  significantly  easier. 

The  problem  with  an  affirmative  answer  which  may  likely  be 

a  "practical"  answer,  to  the  question  posed  above,  is  that  there 

doesn't  seem  to  be  a  theoretical  or  rigorous  analytical  reason 

(to  the  best  of  the  author^'  knowledge)  why  X  /A  must  always  be 

c  c 

approximately  1.  Therefore,  if  on  the  other  hand,  the  answer  to 
our  question  is  negative,  then  the  following  observations  may  be 
made:  (a)  significant  variation  in  k^  is  possible  with  variations 

in  effective  nondimensionalized  coupling  (Ac/Ac).  In  fact  it  can 
be  seen  from  Figure  3  that  if  Ac/Ac  approaches  zero,  k-^  (and  hence 
uj^)  approaches  zero.  (b)  The  values  of  k^  vary  significantly  with 
Ac  for  low  Ac  but  approach  asymptotic  values  for  large  Ac. 

(c)  The  highest  values  of  k^  is  for  isotropic  (metals)  or  quasi¬ 
isotropic  configurations.  (d)  For  large  values  of  A^,  there  ap¬ 
pears  to  be  a  simple  approximate  (hopefully  linear  relationship) 
between  k^  and  ^c/^c  (or  a  measure  of  coupling).,  (e)  For  very 
large  coupling  (Ac/A  -*■()),  k^  approaches  zero,  which  may  provide 
the  ingredient  necessary  for  coupling  between  the  elastic  and  ri¬ 
gid  motions. 


68 


Perhaps  a  number  of  implications  of  some  of  these  observat i  <  mis 
should  be  examined:  Observation  (c)  seems  to  imply  that  the  high¬ 
est  first  frequency  would  correspond  to  isotropic  or  quasi-isotrop- 
ic  configurations  if  (D22/p)  and  £q  are  the  same.  It  is  known  how¬ 
ever  the  metals  have  lower  values  of  (D22/p)  than  composites.  It 
therefore  means  that  quasi-isotropic  or  orthotropic  configurations 
are  desirable  for  such  a  design  goal.  Observation  (e)  would  seem 
to  imply  that  if  a  designer,  interested  in  tailoring  the  wing  fre¬ 
quencies,  arbitrarily  introduces  large  effective  nondimensionalized 
coupling  ( X c/^c~>'0 )  ,  then  (and  hence  m-^)  would  approach  zero. 

This  may  result  in  coupling  between  flexible  and  rigid  motions 
which  may  or  may  not  lead  to  instabilities.  Could  this  have  hap¬ 
pened  in  the  case  of  the  X-29  Forward  Swept  Composite  Wing  Aircraft 
for  which  one  of  the  primary  modes  of  instability  results  from  the 
coupling  between  flexible  and  rigid  body  motions?  In  other  words 
was  "too  much"  coupling  (effective)  inadvertently  built  into  the 
wing  during  the  design  process?  If  that  is  the  case,  is  there  an 
alternative,  equivalent  design  without  any  penalties  (weight  or 
otherwise)  that  could  have  been  explored?  Although  the  answers  to 
these  questions  can,  strictly  speaking,  only  be  possible  after  car¬ 
rying  the  necessary  aeroelastic  analysis  in  which  unsteady  aerody¬ 
namic  forces  are  considered,  it  appears  from  Figure  3  that  a  rough 
idea  of  the  final  picture  may  be  obtained  from  the  natural  fre¬ 
quency  analysis.  Afterall,  it  is  a  common  belief  that  the  phenom¬ 
ena  that  actually  lead  to  aeroelastic  instabilities  are  linked  to 
damping  and  coupling. 
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5.3.9 

Modal  Transformation 

Earlier  in  this  paper  it  was  mentioned  that  previous  studie 
by  other  investigators  have  found  what  appeared  to  be  some  kind 
of  modal  transformations  as  ply  orientation  was  changed  in  a  de- 
sign  process  for  a  composite  wing.  While  variation  in  ply  orien 
tation  may  change  several  directional  stiffness  parameters  for 
the  wing,  the  coupling  stiffness  parameter,  ,  may  be  singled 
out  as  a  significant  design  parameter,  because  it  may  vary  con¬ 
siderably  (orthotropic  conf igurations  have  zero  values  while  it 
may  be  fairly  significant  in  other  configurations) .  Furthermore 
it  should  be  remembered  that  the  main  reason  for  ply  orientation 
variation  is  for  'tailoring',  which  is  believed  to  be  primarily 
tied  to  couplings  (D.^  and  C>2g)  •  The  absence  (or  the  presence) 
of  these  couplings  is  basically  what  differentiates  orthotropic 
configurations  from  anisotropic  configurations.  From  these  ob¬ 
servations,  and  the  fact  that  the  entity  that  ties  the  bending 
and  torsional  equations  is  the  coupling,  it  became  clear  that 
the  role  of  coupling  in  modal  transformations  should  be  signifi¬ 
cant  . 

In  order  to  see  the  role  of  coupling  therefore  in  this  stu¬ 
dy,  the  modal  assumptions  for  the  coupled  problem  were  made  sim¬ 
ilar  to  those  normally  made  for  the  uncoupled  problem  (e.g.,  the 
frequency  was  assumed  to  be  real)  so  as  to  provide  an  opportun¬ 
ity  to  compare,  contrast  and  discern  the  final  results  easily. 
When  this  was  done  and  the  eigen-problem  was  formulated  result¬ 
ing  in  a  complex  transcendental  expression  from  which  the  eigen¬ 
value  are  to  be  extracted,  a  careful  examination  began. 

A  significant  difference  between  the  coupled  and  uncoupled 
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problems  was  that  (as  shown  in  Figure  2)  the  transcendental  ex¬ 
pression  from  which  the  eigenvalues  are  extracted  is  complex  for 
the  coupled  problem  while  it  is  real  for  uncoupled  problems. 

The  complex  nature  of  this  transcendental  expression  basically 
reflects  the  fact  that  the  bending  and  twisting  oscillations  are 
generally  out  of  phase.  Therefore  the  resultant  coupled  frequen¬ 
cy  that  represents  both  the  bending  and  twisting  oscillation  may 
be  viewed  as  some  kind  of  vector  representation  of  the  individual 
contributions.  In  order  to  formulate  some  explanations  for  the 
phenomenon  of  modal  transformation  in  coupled  (conservative)  sys¬ 
tem,  it  may  be  necessary  to  compare  and  contrast  coupled  systems 
with  damped  systems.  Damped  systems,  by  definition  are  noncon¬ 
servative,  i.e.,  the  system  experiences  a  net  loss  or  gain  in  en¬ 
ergy.  It  is  well  known  that  in  a  damped  system  the  transcenden¬ 
tal  expression  for  extracting  the  eigenvalues  is  complex,  again, 
due  to  the  phase  angle  that  exists  between  the  damping  force  and 
the  conservative  forces  in  the  system.  It  is  also  known  that 
some  damping  (desirable  types)  would  tend  to  reduce  the  oscilla¬ 
tion  of  the  system  (the  non-desirable  type  tend  to  make  the  oscil¬ 
lations  diverge) .  Therefore  since  damping  is  linked  to  some  ener¬ 
gy  transfer  which  in  turn  tend  to  lead  to  a  change  in  the  oscilla¬ 
tion  frequencies,  it  was  thought  that  the  complex  nature  which  is 
common  to  the  coupled  and  damped  system  determinants  (transcenden¬ 
tal  functions)  from  which  the  eigenvalues  are  extracted  may  be  a 
similarity  that  may  provide  some  explanations  to  the  modal  trans¬ 
formations  in  coupled  systems.  Using  the  similarity  argument,  the 
coupled  system  which  for  the  present  problem,  is  conservative,  may 


be  viewed  as  having  some  conse 


rvative  inter-modal  energy  trans¬ 
fer  within  the  system  when  the  coupling  is  changed,  resulting 
m  steady  state  changes  or  transformations  of  the  modal  energy 
content  of  a  coupled  mode  compared  to  the  uncoupled  case.  it 
may  be  worthwhile  to  point  out  that  some  results  recently  ob¬ 
tained  at  Purdue  University8  and  communicated  to  the  authors 
seem  to  strongly  support  this  hypothesis. 


5.3.10 

Empirical  Relations 

A  careful  study  of  Figure  3  has  led  the  authors  to  propose 
the  following  closed  form  asymptotic  relationship  that  may  be 
useful  for  some  preliminary  design  consideration: 

k.  =  3.5  (I  /A  )  ,  *  >3-0  (28  ) 

X  O  C  v 

Equation  28  was  derived  from  Figure  4.  Equation  28  as  well  as 
equation  27  show  that  the  first  coupled  frequency  decreases  with 
increasing  coupling,  a  trend  that  seems  to  be  supported  by  new 
results  from  Purdue  University8  and  the  data  from  MIT.3  In  ref¬ 
erence  3  for  example,  the  first  nondimensionalized  frequency  com¬ 
puted  by  Raleigh-Ritz  (in  which  coupling  is  zero)  had  a  value  of 
3.52,  which  is  consistently  higher  than  those  computed  by  finite 
element  method  in  which  coupling  is  finite  (not  equal  to  zero) . 
Equations  27  and  29  which  are  closed  form  (generally  rare  for  an- 
istropic  systems)  should  be  easy  to  use. 


Before  this  discussion  is  concluded,  it  is  probably  necess¬ 
ary  to  explain  why  only  the  results  of  the  first  mode  are  shown 
in  this  paper.  First  of  all  it  should'be  pointed  out  that  some 
second  mode  data  have  been  generated  but  are  still  being  studied 
very  critically  to  understand  the  general  trends.  It  may  also 
be  pointed  out  that  the  extraction  of  the  eigenvalues  is  a  little 
challenging  since  some  care  is  needed  in  ordering  the  roots  of 
the  operator  equations. 


5.3.11 

Concluding  Remarks 

This  paper  has  attempted  to  present  exact  closed  form  solu¬ 
tions  to  the  coupled  bending-torison  vibration  problem  for  a 
simplified  model  of  composite  aircraft  wings  with  war  ping  con¬ 
straint.  Increasing  the  coupling  was  found  to  decrease  the  first 
coupled  frequency.  A  comparison  between  the  coupled  problem  and 
a  sample  damped  problem  led  the  authors  to  propose  some  explana¬ 
tion  to  the  "modal  transformation"  phenomenon  found  by  earlier  in¬ 
vestigators.  Some  simplified  closed  form  expressions  are  provi¬ 
ded  for  the  first  coupled  frequencies  which  may  be  useful  for  fast 
applications . 
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5.4 

1.  Introduction 

The  analytical  study  of  free  vibration  of  a  multi-degree-freedom  system 
presents  itself  in  the  form  of  an  eigenvalue  problem,  in  which  the  eigenvalues 
yield  the  natural  frequencies  and  the  corresponding  eigenvectors  are  the  mode 
shapes.  Together  the  natural  frequencies  and  the  mode  shapes  constitute  the 
"fingerprints"  that  uniquely  characterize  a  dynamic  system.  Moreover,  the 
solu  ion  of  a  forced  vibration  problem  relies  predominantly  upon  the 
technique  based  on  modal  series  expansion  which,  as  implied  by  the  name, 
requires  the  determination  of  eigenvalues  and  modal  vectors.  Indeed,  the 
forced  response  is  obtained  as  a  superposition  of  the  system’s  normal  modes, 
each  one  of  them  being  multiplied  by  a  time-dependent  generalized 
coordinate.  In  the  area  of  dynamic  aeroelasticity,  the  knowledge  of  natural 
modes  is  also  essential.  The  modal  matrix  formed  by  columns  of  eigenvectors 
is  used  as  a  transformation  matrix  to  reduce  the  number  of  degrees  of 
freedom  in  the  flutter  analysis.  By  modifying  the  mode  shapes  of  a  system 
through  aeroelastic  tailoring,  the  flutter  margin  can  be  changed  for  better  or 
for  worse.  It  is  therefore  instructive  to  examine  the  influence  of  the  tailoring 
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parameters  upon  the  vibration  modes.  In  previous  work  tailoring  parameters 
were  suggested  nondimensioned  in  terms  of  combinations  of  the  plate  bending 
stiffnesses.  When  these  parameters  are  varied  within  the  constraint  of  their  j 

applicable  bounds,  the  natural  frequencies  and  mode  shapes  change.  The  j 

magnitude  and  extent  of  these  changes  are  the  subject  of  this  chapter. 

5.4 

2.  Effect  of  Cross  Stiffnesses  on  the  Natural  Frequencies 

To  illustrate  and  to  discern  the  probable  extent  of  the  influence  of 
dimension  parameters  on  plate  natural  frequencies,  a  cantilevered, 
rectangular,  flat  plate  of  differing  aspect  ratios  was  chosen  for  the  study.  The 
geometry  of  the  plate  is  shown  in  Figure  1.  The  deformation  of  this 
planform  is  modelled  by  classical  plate  theory  using  the  Rayleigh-Ritz 
method.  For  this  part  of  the  tailoring  work,  the  assumed  displacement 
function  is  a  polynomial  in  x  and  y  coordinates,  i.c. 

w(x,y)  =  V  V  a^  x1_1  yJ_1  for  m,  n  =  1,2,. ...8  (  L) 

'  =  1  pi 


If  both  the  integers  m  and  n  take  on  the  same  value,  say  4,  then  the 
displacement  function  is  said  to  be  biquartic  (similarly,  bihextal  for 
m  =  n  =  6,  bioctal  for  m  =  n  =  8,  etc...).  For  the  plate  in  Figure  1,  it  has 
one  chordwise  edge  built-in  and  the  rest  of  edges  free.  The  geometric 
boundary  conditions  (displacement  and  slope)  are  enforced  along  y  =  0,  and 
the  displacement  function  is  reduced  to 


m-»  1  n+  1 


w(x,y)  =  V  V  x,_]  yJ_1.  for 

i=l  j=3 


m  =  1,2,.. .,8 

n  =  2,3,. ..,8 


The  advantage  in  using  polynomials  versus  beam  mode  shape  functions  lies  in 
the  relative  ease  of  integration  over  the  trapezoidal  planform  for  which  a 
rectangle  is  a  special  case.  For  the  particular  planform  of  interest,  the  largest 
exponent  order  in  chordwise  displacement  (x-direction)  is  qradratic,  i.e. 
m  =  2.  This  provides  for  a  chordwise  curvature  or  camber  effect  with  a 
parabolic  shape.  The  largest  spanwise  order  (y-d-  •••'tion)  is  quartic,  i.e. 
n  =  4.  This  arrangement  of  terms  is  made  for  two  reasons.  First,  the  plate  is 
supported  along  y  =  0,  thus  restricting  the  freedom  in  chordwise  distorsion. 
Secondly,  the  fewer  number  of  terms  included  in  the  assumed  displacement 
allows  greater  computational  efficiency.  Another  case  study  discussed  later 
will  show  that  no  significant  loss  of  numerical  accuracy  is  entailed  by  using  an 
approximation  of  these  orders. 

Now  applying  Eqn  (2)  to  Ritz  analysis  which  seeks  a  stationary  solution 
to  the  variational  condition  on  the  energy  expression 


6  (V  -  T)  = 


where  V  and  T  are  respectively  the  strain  and  kinetic  energies  of  the  flexural 
laminated  plate,  and  are  written  as  follows: 


Returning  to  our  discussion.  Figures  2  through  9  display  the  non- 
dirnensionalized  frequency  associated  with  the  first  four  modes  plotted  against 
the  cross-stiffness  tailoring  parameters,  V’p  V2  and  The  specimen  considered 
is  the  aforementioned  cantilevered  rectangular  planform  with  varying  aspect 
ratios  (Figure  1).  The  reference  frequency  used  for  the  non- 
dimensionalization,  denoted  as  u.iiref  (where  i  =  1,2, 3, 4),  corresponds  to  that 
which  would  be  computed  with  zero  cross-stiffness  at  each  aspect  ratio  and 
each  ith  mode.  Note  that  the  parameters  d2  and  measure,  respectively, 

’ 1  0  degree  of  camber/twist,  bend/twist  and  camber/bend  interaction  present. 

In  every  figure,  one  of  these  parameters  is  changing,  while  the  other  two  are 

kept  constant.  Care  was  exercised  to  vary  the  tailoring  parameters  within 

their  allowable  ranges.  In  addition  to  the  relation  which  defines  the 

parameter’s  upper  and  lower  bounds,  another  constraint  was  also  satisfie< 

Noncompliance  with  these  constraints  would  not  cause  computational  troubles 

in  most  cases,  but  it  would  deprive  the  analysis  of  any  physical 

meaningfulness. 

Observing  the  figures,  one  can  discern  the  general  trend  that,  for  small 
values  of  the  tailoring  parameters  (e.g.  below  0.3  ~~  0.4),  stiffness  cross- 
coupling  has  little  influence  on  the  natural  frequencies.  As  the  parameters 
approach  their  bounds,  then  the  dependence  of  the  frequencies  on  the 
coupling  is  extremely  evident  for  most  cases.  Frequencies  generally  decline 
significantly  as  the  couplings  accentuate  their  effects,  although  there  are 
several  instances  showing  opposite  behavior  to  this  trend.  In  the  cases  where 
frequencies  do  increase  with  stronger  couplings,  the  magnitude  of  the  increase 
is  not  as  substantial  as  that  of  the  decrease  with  larger  couplings.  Compared 
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to  previous  investigations  by  Weisshaar,  frequencies  of  a  beam  model  were 
found  to  be  relatively  unaffected  in  small  couplings,  but  decline  when  t/' 
attains  high  values.  Furthermore,  in  Figures  2  and  3,  one  does  not 

observe  any  symmetry  of  the  frequency  curves  with  respect  to  i/!j  =  0. 
However  slightly  distorted,  this  symmetry  begins  to  take  shape  when  tv,  and 
are  variables,  as  depicted  in  Figures  4  through  7.  When  one  of  the  elastic 
couplings  is  removed,  such  as  camber/twist  interaction  in  Figures  8  and  9, 
then  the  frequency  curve  symmetry  with  respect  to  i;2  =  0  is  perfect.  Finally, 
as  seen  in  the  pattern  of  variations,  no  1  ^rvation  can  be  drawn  on  the 
influence  of  the  aspect  ratios. 

5.4.3 

3.  Effect  of  Cross  Stiffnesses  on  the  Mode  Shapes 

Of  key  interest  to  an  aeroelastician  is  the  effect  of  structural 
deformations  of  a  lifting  surface  upon  the  aerodynamic  loads.  In  the  dynamic 
environment,  the  surface  deformations  associated  with  the  mode  shapes  exert 
a  strong  bearing  on  the  flutter  stability.  Hence,  tailoring  for  desirable  modal 
shapes  is  a  passive  flutter  prevention  technique.  The  most  significant 
alteration  of  a  mode  shape  is  the  repositioning  of  nodal  line  contours.  A 
nodal  line  is  a  locus  of  zero  displacement  present  on  the  structure  when  that 
structure  vibrates  at  a  particular  natural  frequency.  Biot  and  Arnold 
reported  computations  wherein  zero  airspeed  flutter  occurs  for  a  two- 
dimensional  airfoil  when  a  node  (point  in  this  case)  is  located  at  the  3/4  chord 
position  aft  of  the  leading  edge.  They  also  stated  that  in  a  vortex-free  flow, 
i.e.  flow  without  circulation,  flutter  occurs  at  all  airspeeds  when  a  node 
coincides  with  the  3/4  chord  position.  In  this  latter  case,  frequency  merging 
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takes  place,  and  the  ratio  of  pure  bending  to  pure  torsion  frequencies  is  nearly 
unity.  The  authors  concluded  that,  in  designing  an  airfoil,  it  is  recommended 
to  place  the  nodal  lines  as  far  away  as  possible  from  the  3/4  chord  "danger 
zone".  This  is  an  illustration  of  the  importance  of  mode  shapes  to  flutter 
occurrence. 

The  potential  for  modifying  mode  shapes  through  stiffness  cross-coupling 
is  illustrated  in  Figures  10  through  12.  In  these  figures,  the  modes  are 
represented  in  terms  of  nodal  lines.  Using  the  same  specimen,  nodal  fines  for 
the  first  six  mo  .  of  the  fixed-aspect-ratio  planform  are  reproduced  for  a 
number  of  values  of  vt  and  U2.  These  nodal  lines  are  generated  from  the 
Rayleigh-Ritz  analysis.  In  each  of  the  figures,  two  of  the  three  cross-coupling 
parameters  il\,  r2  and  are  held  constant,  while  the  third  ranges  within  the 
constraints  imposed.  Note  that  the  spanwise  length  of  the  rectangular 
planform  has  been  scaled  down,  so  that  displacement  data  in  one  mode  shape 
are  displayed  on  one  computer  output  page.  Two  comments  are  in  order.  The 
different  modes  are  not  symmetric  with  respect  to  il\  —  0.0  in  Figure  10. 
However,  in  the  following  two  figures  where  il'2  is  varying,  modal  symmetry  is 
clearly  discernible.  As  a  matter  of  fact,  this  symmetry  appears  perfect  in 
Figure  12  where,  again,  the  camber-twist  modes  are  decoupled.  The  above 
situation  is  consistent  -  ith  that  observed  for  the  frequency  curves.  Based  on 
the  observations  made,  one  can  therefore  infer  that  the  camber/twist 
interaction  can  behave  very  differently  when  the  sign  of  the  elastic  coupling 
changes.  This  is  in  sharp  contrast  to  the  other  two  types  of  couplings. 

As  the  second  comment,  a  reordering  of  the  modes  occurs  between,  for 
instance,  the  second  and  third  modes  in  Figure  10.  The  reordering  occurs 


when  C'i  changes  its  value  from  0.6  to  0.8.  At  i/'j  =  0.6,  the  second  mode  is 
predominently  a  second  bending  mode,  while  the  third  mode  is  predominently 
a  first  torsion  mode.  At  4:\  =  0.8,  the  second  mode  becomes  predominently  a 
first  torsion  mode,  whiie  the  third  mode  predominently  a  second  bending 
mode.  This  change  in  modal  identification  is  sometimes  referred  to  as  a 
modal  exchange  or  modal  transfer.  The  tendency  for  the  two  modes  in 
question  to  interchange  their  order  can  also  be  seen  in  the  extreme  negative 
range  of  t in  Figure  10.  In  Figure  11  and  12,  modal  exchange  involves 
the  fourth  and  fifth  modes  as  increases.  In  all  the  exam cited,  modal 
exchange  occurs  as  the  result  of  changes  in  one  of  the  cross-stiffnesses  or 
t'2).  This  is  important  to  dynamic  behavior  because,  while  the  modal 
frequency  is  relatively  unaffected  by  changes  in  these  parameters,  the  modal 
forces  (also  called  generalized  forces),  that  arise  as  a  result  of  deformation  in 
that  mode,  may  be  affected  significantly. 

The  need  now  arises  to  analyze  and  to  explain  the  phenomenon  of  modal 
transfer  in  a  more  quantitative  manner.  Oyibo  in  a  study  funded  by 
AFOSR  Contract  F49620-85-C-0090 ,  first  proposed  that  this  phe¬ 
nomenon  can  be  explained  in  terms  of  inter-modal  energy  trans¬ 
fer.  In  the  next  sections  an  independent  approach  shall  be 
used  to  try  to  verify  Oyibo’ s  theoretical  explanation. 

5.4.4 

4.  Modal  Content  Analysis 

For  a  vibrating  anisotropic  plate,  the  mode  shapes  are  com¬ 
posed  of  a  number  of  intricately  coupled  fundamental  shapes. 

As  a  result  of  the  presence  of  cross-coupling,  the  distinct  ca¬ 
tegorization  of  coupled  modes  as  pure  bending  or  pure  torsion 
modes  is  not  possible.  A  previous  work  discusses  the  idea  of 
expression  of  mode  shapes  for  a  c ross-coupled  structure  in 
terms  of  a  set  of  modes  generated  from  a  baseline  configuration. 
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This  decomposition  may¬ 
be  accomplished  in  terms  of  so-called  modal  participation  factors  as  explained 
in  the  following  discussion. 

Consider  the  equation  of  motion  for  free  vibration  of  a  special  plate, 
called  the  baseline  structure.  This  baseline  structure  has  no  stiffness  cross¬ 
coupling  between  the  bending  and  torsion  modes  (i.e.  f.'j  =  (.'2  =  0).  The 
equation  is  written  as: 

-  ~o2  (Mol  M  +  (Koi  {’/I  -  (0)  (7) 

where  |M0]  and  [K0j  are,  respecth  c  i.  ,  the  mass  and  stiffness  matrices  of  the 
elastically  uncoupled  system.  The  solution  to  the  above  eigenvalue  problem 
provides  the  modal  vectors  {r/1)}  of  the  elastically  uncoupled  system.  The  full 
modal  matrix  Jr;],  a  square  matrix,  can  be  constructed  by  assembling  columns 
of  eigenvectors  {7/^}  in  ascending  order  of  the  natural  frequencies.  Consider 
now  a  plate  with  nonzero  cross-couplings.  Let  [K]  be  the  portion  of  the 
stiffness  matrix  due  to  the  effects  of  rotating  the  ply  angle  to  create  ‘iffness 
cross-coupling.  This  matrix  may  be  added  to  the  uncoupled  portion  (K0]  of 
the  baseline  structure.  As  a  result,  we  encounter  a  new  eigenvalue  problem, 
for  which  V’i  and  t/:2  are  nonzero,  as  follows: 

-  J  |M0|  {?}  +  IK„+K|  {0  =  {0}  (8) 

The  eigenvalues  and  eigenvectors  of  the  free-vibration  problem,  Eqn  (5.7),  are 
different  than  those  found  in  Eqn  (8)  because  of  the  modified  stiffnesses.  A 
coordinate  transformation  using  [rj]  may  be  performed  to  obtain  the  following 
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approximation  for  the  new  modal  vectors  {t}  of  the  elastically  coupled 
system: 

{chh  =  [,/j  (phh  (9) 

The  ith  element  of  {P^},  denoted  as  P  ,  is  the  contribution  of  the  ith 
uncoupled  mode  to  the  jth  coupled  mode  {£^}.  Thus  we  see  that  the  jth 

modal  vector  of  the  coupled  structure,  {£^},  is  a  linear  combination  of  the 
modal  vectors  of  the  uncoupled  baseline  structure,  the  coefficients  in  the 
combination  being  the  elements  of  {P^}.  This  vector  is  called  the  moda 
content  vector.  Now  inserting  Eqn(9)  into  Eqn  (8),  premultiplying 
throughout  by  [r/jT,  and  normalizing  the  modal  amplitude  so  that  the 
following  two  relationships  hold: 

lnlT  |M„|  In]  =  !  i  ] 

|r,]T  [KJ  I.,!  =  UM 

where  [  I  j  is  the  identity  matrix  and  [Tu02.J  a  diagonal  matrix  whose  nonzero 
elements  are  the  natural  frequencies  of  the  baseline  structure.  We  now  find 
that  Eqn  (8)  becomes 

-4*  {P,J)}  +  (V-l  +  |nlT  |K|  |n|  {P,M>  =  {o}  (12) 


If  we  further  define 


then  we  can  rewrite  Eqn  (12)  as 


{P(jl!  +  I'-bM  {P1’1!  +  |K|  =  {0}  (  H) 

To  determine  the  elements  of  the  modal  content  vector  \P^}  associated  with 
the  jth  coupled  mode,  we  begin  first  by  solving  Eqn  (7)  ,  finding  u.02  and  {r/} 
for  all  uncoupled  modes.  Next  tv'2  of  jth  coupled  mode  is  determined  from  the 
eigenvalue  problem,  Eqn  (8)  .  Finally  the  vecti  ^P^lj,  which  is  itself  an 
eigenvector  of  Eqn  (12)  ,  can  be  calculated  from  the  system  of  homogeneous 
equations,  Eqn  (12)  if  one  of  its  elements  is  arbitrarily  rhosen.  There  is, 
however,  a  more  convenient  way  to  obtain  {P^}.  Solving  for  all  the  modal 
vectors  {r?^}  (where  i  =  1,2, ...,n)  and  the  jth  mode  of  Eqns  (7)  and 

(5.8)  respectively,  premultiplying  Eqn  (9)  by  [r)]T  [M0]  and  normalizing  in 
accordance  with  Eqn  (10)  we  get 

{P(J)}  =  [l]T  [M0]  {$<’>}  (  15) 

If  we  normalize  {P^}  by  dividing  its  elements  by  the  largest  absolute  value  in 
{P^},  then  we  get  the  relative  contribution  of  each  mode  (r/1)}  to  make  up 
the  content  of  each  {£^}.  Using  the  above  technique  to  determine  modal 
content,  we  can  pinpoint  those  modes  that  may  be  important  to  the  flutter 
analysis. 


The  modal  content  vector  {P^}  can  also  be  applied  to  the  energy 
approach  of  the  vibration  problem.  Consider  once  more  the  elastically  coupled 


free  vibration  equation: 

-  iXi  if*'1)  +  |K0+K|  if*'1)  =  {0}  (8) 

Substituting  Eqn  (9)  -  the  relationship  between  the  uncoupled  and  coupled 

eigenvectors  -  into  Eqn  (8)  and  proceeding  as  before,  we  obtain  Eqn  (12)  . 
Then,  premultiplying  Eqn  (12)  throughout  by  the  row  vector  tP^\,  we  get 

-,-j2  P!j),  {P^}  +  [P(j),  ['a02J  {PM}  +  tPW,  [K]  {P(J)}  =  {0}  (16) 

At  this  point,  we  can  form  the  Rayleigh’s  quotient  from  Eqn  1 16)  as  follows: 

^ ,  _  .p(il.  iV-l  {p|j)?  +  iP(j)i  H  {p*jI>  U7) 

"J  tP(j),  {P(j)} 

By  adopting  the  notation 


{PM}  = 


\ 

PjU) 

pu 

p.w 

*- 

• 

•  =  ■ 

p2j 

p  (j) 

A  n 

pBJ 

; 

we  can  rewrite  Eqn  (17)  in  summation  form.  i.e. 


(18) 


an 


£  -o,2  p,/  +  £  £  Klk  p,j  pkJ 

i-l  i-l  k-i 


y  P  1 

Lj  1  1J 
1  =  1 


(19) 


For  the  jth  mode  of  the  elastically  coupled  system,  the  numerator  in  Eqn 
(19)  represents  the  strain  energy.  The  first  term  can  be  viewed  as  the 
primary  strain  energy  of  the  baseline  structure,  while  the  second  term  is 
interpreted  as  the  strain  energy  due  to  stiffness  cross-coupling.  The 
denominator  represents  a  measure  of  the  kinetic  energy. 

In  the  previous  section,  we  observed  the  phenomenon  of  modal  transfer. 

To  examine  this  phenomenon  from  an  energy  viewpoint,  we  focus  our 

n 

attention  on  the  kinetic  energy  expression  P,j  .  First  we  define  each 

i-l 

element  P^2  of  [P2!  matrix  as  the  contribution  of  the  kinetic  energy  of  the  ith 
uncoupled  mode  of  the  baseline  structure  to  the  total  kinetic  energy  of  the 
mode  of  the  coupled  structure.  To  facilitate  further  the  numerical 
interpretation  of  P,,2,  each  ith  element  of  the  jth  column  of  [P2]  is  divided  by 
the  sum  of  the  elements  in  that  column.  As  a  result,  we  get  the  sum 

n 

V  Ptj2  =  l  for  the  jl  column  and  obtain  a  normalized  energy  contribution  of 

i  =  i 

each  ith  uncoupled  mode  to  that  of  the  jtk  coupled  mode. 

As  an  illustrative  example,  a  cantilevered  laminate  of  rectangular 
pianform  with  the  stacking  sequence  [#2/0js  was  chosen.  The  plies  are  of  equal 
thickness;  the  ply  orientation  6  is  measured  from  the  longitudinal  axis  and 
takes  on  values  from  0  to  90  degrees.  The  baseline  structure  is  defined  as 
02/0ls;  it  is  orthotropic  (D16  =  D26  =  0)  and  elastically  uncoupled.  There  are 
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6~  21  5  - 


0.99  0.01  0 

0.01  0.49  0.49 
0  0.49  0.48 

0  0.01  0.02 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 


0  0  0 
0  0  0 
0.03  0  0 

0.83  0.11  0 

0.09  0.68  0.20 
0.04  0.02  0.28 
0.01  0.01  0.14 
0  0.17  0.35 

0  0.01  0.01 


0  0  0 
0  0  0 
0  0  0 
0.02  0.01  0 
0  0.02  0 
0.62  0.04  0 

0.23  0.61  0 

0.12  0.29  0.06 
0  0.04  0.94 


20) 


As  f>  increases  further,  the  bending  kinetic  energy  begins  to  predominate, 
while  the  torsional  kinetic  energy  gradually  declines.  Finally,  near  8  —  65  °  , 
mode  2  becu:,t.~  a  nearly  pure  flexural  mode  (P3-,2  =  l).  Similar  comments  can 
be  made  with  respect  to  mode  3  in  Figure  14  where  it  is  the  second  pure 
bending  mode  at  8  =  0  0  (P332  =  1)  and  transmutes  into  the  first  pure  torsion 
mode  at  8  =  65  0  (P232  =  l)-  Thus,  through  an  evolutionary  process  which 
consists  of  changing  gradually  the  stiffness  characteristics,  we  observe  a 
reciprocal  switch  or  exchange  of  modal  character,  and  in  the  context  of 
kinetic  energy,  this  modal  exchange  is  rendered  possible  by  an  exchange  of 
modal  kinetic  energies. 

For  higher  modes,  say  and  {£^},  the  deformation  pattern  is  more 

complex.  As  a  result,  there  are  a  larger  number  of  uncoupled  modes 
participating  in  the  motion  of  the  coupled  system.  Hence,  modal  content  is 
more  distributed,  as  is  depicted  in  Figures  16  and  17.  A  further  remark  is 
in  order  here.  The  fact  that  the  coupled  modes  {^}  and  {^}  involve  only 
the  interaction  of  first  torsion  and  second  bending  is  well  demonstrated  by  the 
equalities  P322  =  P232  and  P222  =  P332  in  Figures  13  and  14.  Such 
equalities  do  not  exist  in  Figures  16  and  17.  For  {£^}i  there  are  three 
modes  of  the  baseline  structure  that  are  coupled  together  when  6  #  O'  ,  while 
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two  reasons  for  selecting  this  particular  laminate:  (a)  the  layup  is  unbalanced, 
hence  it  has  a  greater  inherent  coupling  of  deformations  as  compared  to  the 
balanced  counterpart;  (b)  the  dynamic  behavior  of  this  particular  laminate 
has  been  documented  previously  by  Jensen,  et  al.  [MIT].  This  permits 
correlation  of  our  study  to  established  results.  Another  advantage  to  using 
this  unbalanced  laminate  is  that  modes  2  and  3  involve  only  an  interplay 
between  the  first  torsion  and  second  bending  behavior  of  the  baseline 
structure.  The  number  of  terms  considered  in  Eqn(2)  were  again  set  to 
m  =  2  and  n  =  4,  resulting  in  a  nine-c  °e-of-freedom  system  (see  Eqn  (2)  ). 
These  orders  of  the  exponent  yield  numerical  results  in  satisfactory  agreement 
with  those  obtained  experimentally  by  Jensen;  this  is  demonstrated  in  Figure 

to  be  commented  shortly.  Now  referring  to  Figure  13  and  starting  at 
ti  =  0°  ,  mode  2  is  observed  as  having  all  the  kinetic  energy  of  the  pure  first 
torsion  mode  (P222  =  1  at  8  =  0 0 );  there  is  no  contribution  of  the  kinetic 
energy  from  the  second  bending  mode  (P322  =  0  at  6  *=  0  ° ).  The  nodal  line 
patterns  at  6  =  0  0  are  sketched  in  Figure  15  along  with  those  corresponding 
to  other  ply  angles. 

In  Figure  13,  it  is  noted  that,  as  6  is  increased,  mode  2  becomes  more 
coupled;  there  occurs  a  gradual  decrease  of  the  torsional  modal  kinetic  energy 
P2n2  with  the  corresponding  increase  of  the  bending  modal  kinetic  energy  P322. 
At  about  6  =  21.5*,  the  modal  energy  curves  cross  each  other,  and  mode  2 
has  49%  of  first  torsional  energy  and  another  49%  of  second  bending  energy, 
the  remainder  being  contributed  by  other  fundamental  modes.  For  purpose  of 
illustration  of  the  modal  energy  distribution,  the  normalized  kinetic  energy 
matrix  for  0  =  21.5*  is  shown  as  follows: 
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for  {i^},  there  are  four  modes.  Because  of  the  same  number  and  type  of 
fundamental  modes  participating  in  the  motion,  the  crossings  of  the  energy 
curves  in  Figures  13  and  14  occur  at  the  same  angle  8  and  indicate  the 
occurrence  of  a  modal  transfer.  If  the  modal  participation  is  different,  as  it  is 
the  case  in  Figures  16  and  17,  then  the  crossings  of  the  energy  curves  do 
not  take  place  at  the  same  angle  8,  however  they  still  suggest  the  occurrence 
of  modal  transfer  somewhere  in  the  vicinity  of  the  crossing  angles  8. 

We  now  turn  attention  to  the  behavior  of  the  natural  frequencies  as  the 
plies  rotate.  Of  interest  in  Figure  18  is  the  tendency  the  frequencies  of 
the  second  and  third  modes  to  approach  each  other  in  the  neighborhood  of 
fl  =  21.5‘.  However  the  frequency  curves  do  not  cross  each  other  and  are 
observed  never  to  do  so.  This  non-crossing  behavior  of  the  frequency  curves  is 
even  more  dramatically  illustrated  by  Ritchie,  et  al.  and  Weisshaar,  et  al. 
In  particular,  Ritchie  mentions  a  so-called  "non-crossing  rule"  in 

conjunction  with  the  frequency  curves.  Ritchie  states  that  the  non¬ 
crossing  rule  forbids  th^  tuning  of  coupled  oscillators,  whether  mechanical  or 
electrical,  to  the  same  frequency.  Although  efforts  were  made  by  the  present 
author  to  find  all  the  stipulations,  if  any,  associated  with  this  rule,  the  search 
has  revealed  only  that  the  term  "non-crossing  rule"  is  well  established  in 
physics.  Specifically,  it  refers  to  the  non-crossing  of  electrical  potential  versus 
internuclear  distance  in  the  context  of  two  electronic  states .  However  we 
are  unable  so  far  to  relate  this  rule  with  either  mechanical  or  electrical 
oscillators  found  in  the  literature.  It  appears  that  Ritchie  used  the 
terminology  in  a  very  loose  manner  and  that,  in  lieu  of  a  rule,  he  really  meant 
the  "non-crossing  behavior". 
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On  the  other  hand,  if  we  compute  the  natural  frequencies  of  the 
structure  versus  the  fiber  orientation  (see  Figure  19),  but  only  take  into 
account  values  of  Dn,  D12,  D22  and  D66  and  disregard  D)6  and  D26,  then  the 
frequency  curves  of  second  and  third  modes  do  cross  at  6  =  21.5”  .  Although 
such  a  system  is  uncoupled,  it  is  not  the  baseline  structure  defined  earlier. 
Instead,  the  present  configuration  refers  to  the  structure  that  is  always 
elastically  uncoupled  by  intentionally  setting  D16  =  D26  =  0,  or  tj  =  v2  =  0 
for  all  fi's.  For  the  baseline  structure  we  have  D16  =  D26  =  0  only  for  6  =  0  0 . 

In  Figure  19,  the  crossing  c  °  uncoupled  frequency  curves  indicates  that  a 
modal  transfer  or  reordering  takes  place.  The  same  can  be  said  with  respect 
to  Figure  20  where  the  double  crossings  between  the  fourth  and  fifth  modes 
suggest  two  modal  transfers.  An  examination  of  the  nodal  patterns  confirms 
the  above  remark. 

In  summary,  we  can  draw  several  conclusions,  based  on  the  above 
discussion: 

i.  Modal  transfer  of  an  elastically-coupled  laminated  plate  can  be  viewed 

and  charted  by  examining  baseline  modal  kinetic  energies;  as  first  pro¬ 
posed  by  Oyibo  f AFOSR  Contract  F49620-85-C-0090*]  . 

ii.  The  frequency  curves  of  two  modes  plotted  against  ply  angle  can 
approach  each  other,  but  appear  never  to  cross  each  other  as  evidenced 
by  a  number  of  plots  of  natural  frequency  versus  ply  angle; 

iii.  The  ply  angle  at  which  modal  transfer  occurs  can  be  deduced  by  the 
crossing  of  the  plots  of  frequency  versus  ply  angle  of  the  elastically 
uncoupled  system  (D16  =  D26  =  0  for  all  6's). 
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Graphite/Epoxy 

Laminate:  [90/135/45/180Js 

=  79.04  lb-in  Density:  p  =  0.055  lb/in3 

$22  =  lb-in  Ply  thickness:  t  =  0.0052  in 

066  =  14.14  lb-in 

D12  =  12.64  lb-in 

D i &  =  4.78  lb- in 

=  4.78  lb-in 


Figure  1.  Geometry  and  Material  Properties  of  a 
Laminated  Plate 


96 


Non-Dimens ional  ized  Frequency,  ,j'(;/|J,2ref  Non-Dimensional  ized  Frequency, 


Camber/Twist  Measure,  ^ 

Figure  2.  Frequency  Variation  as  a  Function  of  Camber/Twist 
Tailoring  Parameter,  in  First  and  Second  Modes 
with  Different  Aspect  Ratios  (AR) 
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Figure  3.  Frequency  Variation  as  a  Function  of  Camber/Twist 
Tailoring  Parameter,  in  Third  and  Fourth  Modes 
with  Different  Aspect  Ratios  (AR) 
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Bend/Twist  Measure,  ^ 

Figure  5.  Frequency  Variation  as  a  Function  of  Bend/Twist 

Tailoring  Parameter,  02,  in  Third  and  Fourth  Modes 
with  Different  Aspect  Ratios  (AR) 
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Benci/Camber  Measure,  y 

Figure  .  6.  Frequency  Variation  as  a  Function  of  Bend/Camber 
Tailoring  Parameter,  '7,  in  First  and  Second  Modes 
with  Different  Aspect  Ratios  (AR) 
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Figure  10.  Effect  of  Tailoring  Parameter  »/',  on  Nodal  Patterns  and 
Natural  Frequencies.  Spanwise  Length  Is  Scaled  Down. 
(AR  =  3,  </',  --  0.273,  1  -  0.303) 
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Figure  11.  Effect  of  Tailoring  Parameter  i/'2  on  Nodal  Patterns  and 
Natural  Frequencies.  Spanwise  Length  Is  Scaled  Down. 
AR  =  3,  i/’i  =■■  0.143,  l  =  0.305) 
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Effect  of  Tailoring  Parameter  V’2  on  Patterns  and 

Natural  Frequencies.  Spanwise  Length  is  Scaled  Down. 
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Figure  13.  Change  of  the  Measure  of  Kinetic 

Second  Coupled  Mode  with  Ply  Angle 


e  (degrees) 

Figure  14.  Change  of  the  Measure  of  Kinetic  Energy  of 
Third  Coupled  Mode  Vith  Ply  Angle 
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Figure  15.  Nodal  Patterns  and  Frequencies  of  the  First  Three 
Modes  Changing  with  Ply  Angle  for  [02/°ls 
Laminate.  Spanwise  Length  Is  Scaled  Down. 
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Fourth  Coupled  Mode  -with  Ply  Angle 


Figure  17.  Change  of  the  Measure  of  Kinetic  Energy  of 
Fifth  Coupled  Mode  with  Ply  Angle 
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Figure  18.  Variation  of  the  Natural  Frequencies  of  the 
First  Three  Modes  with  Ply  Orientation  in 
the  Elastically  Coupled  Laminate 
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Figure  19.  Variation  of  the  Natural  Frequencies  of  the 
First  Three  Modes  with  Ply  Orientation  in 
the  Elastically  Uncoupled  Laminate 
(^1  =  ^=0  for  all  6) 
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'  physical  and' affine  space 
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=  cliordwise  integrals 

“  differential  aerodynamic  pressure  dis¬ 
tributions  in  physical  and  affine  space 
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pararae ters 

»  affine  space  running  aerodynamic  lift 
and  moment,  respectively 
<=■  virtual  work  expressions  in  physical,  and  ■ 
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=>  elastic  constants 

a  affine  space  material,  and  air  density, 
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=  displacement 
=  wing  box  depth 

=  affine  space  bending  and  torsional  dis¬ 
placement,  respectively 
=  affine  space  half-chord  and  chord,  res¬ 
pectively 

=  affine  space  half-span  for  the  wing 
=  parameter  that  measures  the  location  of 
the  elastic  axis  relative  to  mid-chord 
=  affine  space  mass  per  unit  span 
=  vibration  frequency 


^  ^  Introduction 

The  performance  of  modern  supermaneuverable  aircraft  can  be  made 
to  benefit  a  great  deal  from,  significant  advances  in  materials 
technology  and  the  availability  of  more  accurate  aerodynamic  pre¬ 
diction  capabilities.  Supermaneuverability  as  a  design  goal  in¬ 
variably  calls  for  an  optimization  of  the  design  parameters.  Op¬ 
timization  may  be  partially  accomplished  for  example,  by  using 
composite  materials  to  minimize  weight.  Indeed,  it  has  been  known 
that  these  composite  materials  can  be  tailored  to  resolve  the  dy¬ 
namic  or  static  instability  problems  of  these  types  of  aircraft. 
The  concept  is  referred  to  as  aeroelastic  tailoring. 


While  aeroelastic  tailoring  has  tremendous  advantages  in  the  design 
of  an  aircraft,  the  analysis  which  provides  the  basis  for  the  aero¬ 
elastic  tailoring  itself  is  generally  very  involved.  This  is  ra¬ 
ther  unfortunate  since  a  good  fundamental  physical  insight  of  the 
tailoring  mechanism  is  required  for  accurate  and  reliable  results; 
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In  this  paper  an  attempt  is  made  to  look  at  some  dynamic  theories 
that  can  be  used  to  understand  the  aeroelastic  tailoring  mechan¬ 
ism.  Specifically,  the  accuracy  of  the  St.  Venant  torsion  theory 
which  is  relatively  simple  and  frequently  used  in  aeroelastic 
analysis  is  examined  with  particular  reference  to  the  effects  of 
the  wing's  aspect  ratio  as  well  as  other  design  parameters. 

12  3 

Although  earlier  studies  '  '  have  indicated  that  the  St.  Venant's 
torsion  theory  is  reasonably  accurate  except  for  aircraft  wings 
with  fairly  low  aspect  ratios,  the  theory  supporting  that  conclu¬ 
sion  was  based  on  the  assumption  that  the  wing  is  constructed  of 
isotropic  materials.  Basically,  the  St.  Venant' s  torsion  theory 
assumes  that  the  rate  of  change  of  the  wing's  twist  angle  with 
respect  to  the  spanwise  axis  is  constant  (for ; constant  stiffness 
and  torque)  .  This  assumption  is  hardly  accurate  particularly  for 
modern  aircraft  construction  in  which  different  construction  mate¬ 
rials  are  employed  and  the  aerodynamic  loads  vary  significantly 
along  the  wing's  span.  However,  References  1-3  have  shown  that 
(in  spite  of  such  an  inaccurate  assumption)  the  main  parameter 
that  determines  the  accuracy  of  the  St.  Venant 's  theory  is  the 
wing's  aspect  ratio.  Thus,  it  was  determined  that  the  theory  is 
fairly  accurate  for  moderate  to  high  aspect  ratio  wings  construc- 
ted  of  isotropic  materials.  In  recent  studies  '  however,  it  has 
been  shown  that  for  wings  constructed  of  orthotropic  composite 
materials,  the  conclusions  of  References  1-3  need  to  be  modified. 
Rather  than  using  the  geometric  aspect  ratio  of  the  wing  to  deter¬ 
mine  the  accuracy  of  St.  Venant 's  twist  theory,  it  was  suggested 
that  a  generic  stiffness  ratio,  as  well  as  an  effective  aspect 
ratio  which  considers  the  wing's  geometry  and  the  ratio  of  the 
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principal  directional  stiffness,  should  be  considered  in  estab¬ 
lishing  the  accuracy  of  St.  Venant's  theory. 

The  present  paper  is  basically  an  extension  of  the  studies  that 
were  begun  in  References  4  and  5.  In  this  study  the  task  was  to 
examine  the  role  of  coupling  (elastic  coupling)  on  the  accuracy 
of  St.  Venant's  theory  applied  to  static  problems.  It  was  dis¬ 
covered  that  coupling  plays  a  very  significant  role  on  the  ac¬ 
curacy  of  St.  Venant's  twist  theory. 

5.5.3 

Formulation 

Consider  an  aircraft  wing  fabricated  of  composite  materials  and 
mathematically  idealized  as  a  canti-levered  plate  subjected  to 
forces  and  moments.  It  can  be  shown  that  the  equations  of  motion 
for  such  a  model  can  be  described  as  follows. 

alho1V  +  a2ao1V  +  a5ao111+  poalho  + 


p  a0a  =  L 
o  2  o  o 


a  h  *v-ach  aaa+  a-,a  1V-  a. a11  +  p  a^a  + 
2  o  5  o  3  o  4  o  '■>  <  ° 
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)  Cl  t  U 

o  3  o 
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o  2  o  o 
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D. .  are  the  elastic  constants,  p,  is  the  material  density,  Ap  is 
13 

the  differential  pressure  distribution,  w  is  the  displacement,  t 
is  the  time,  and  h  is  the  wing  box  depth. 


5.5.4 

Evolution  of  Warping  Parameters 

The  evolution  of  the  warping  parameter  with  which  to  study  the 
aeroelastic  warping  constraint  phenomenon  for  wings  fabricated 
of  composite  materials  is  a  process  that  depends  on  the  sophisti¬ 
cation  of  the  wing's  mathematical  model;  whether  coupling  ef¬ 
fects  are  included,  whether  the  wing's  chordwise  curvatures  are 
included  and  so  on.  Therefore,  any  warping  parameter  is  as  good 
as  the  corresponding  wing's  displacements  assumptions.  However, 
the  virtual  work  equation  makes  it  possible  for  the  analyst  to 
determine  its  effective  independent  variables  even  before  the  dis- 
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placement  assumptions  are  made.  By  non-dimensionalizing  the 
spanwise  space  variable  in  Equation  (1) ,  depending  on  whether  one 
is  interested  in  the  static,  dynamic,  coupled  or  uncoupled  dis¬ 
placements,  one  of  the  following  warping  parameters  may  be  useful. 


where 

Do  =  D*(1_c)  (7) 

U0/c0)  is  defined  as  the  wing's  effective  aspect  ratio  and  D* 
and  L  are  the  generalized  stiffness  and  coupling  ratios  respective¬ 
ly  (defined  in  earlier  work  such  as  References  5  and  6) . 

Equations  (5)  and  (6)  represent  the  appropriate  warping  parameter 
for  dynamic  deformation,  static  displacement  with  elastic  cross¬ 
coupling  . 

It  was  discovered  in  this  study  that  evolving  the  warping  param¬ 
eter  in  a  manner  shown  in  Equation  (1)  thru  (3) ,  should  enable  one 
to  effectively  investigate  the  effects  of  warping  on  the  composite 

wing's  dynamics  (or  the  accuracy  of  St.  Venant's  theory.)  From  the 

* 

lamination  theory  for  composites  it  is  known  that  while  D0  and 
(i,  /cQ)  are  always  positive.  Lean  be  positive  or  negative. 

However,  from  Equations  (1)  and  (3) ,  it  is  clear  that  whether  a 
composite  wing  has  positive  or  negative  coupling,  the  warping  ef¬ 
fect  (in  terms  of  X  )  is  unchanged. 


5.5.5 


Computations 


By  using  the  evolved  warping  parameters  defined  in  Equations  7  and 
8  and  appropriate  boundary  conditions,  the  boundary  value  problems 
associated  with  Equation  (4)  are  solved  in  a  closed-form  manner  to 
determine  the  wing's  static  twist. 


The  wing  loading  conditions  considered  in  this  analysis  are  as 
follows:  (a)  steady  state  distributed  twist  loads  and  (b)  steady 

state  concentrated  twist  oads. 


5.5.6 


a.  Steady  distributed  twist  loads 

For  a  wing  with  a  constant  uniformly  distributed  spanwise 
twisting  moment,  f  resulting  from  a  steady  state  coupled 
bending-torsion  displacements,  the  exact  closed  form  solutions 
for  the  mode  shape  aQ,  satisfying  the  appropriate  boundary  con¬ 
ditions  is  given  by 
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6  Vo 
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sinh4  A  y 
cJ  o 
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— —  (tanh4A  +  — — - — — )  (cosh4A  y  -1)  I 

4 A  C  4 A  cosh4 A  C  ° 

c  c  c  J 

Equation  8  is  therefore  a  closed  form  coupled  twist  distribution 

for  a  composite  wing  with  the  warping  effects  accounted  for. 

where 


y  =  y  /l  (9) 

o  1  o  o 

When  equation  8  is  evaluated  at  the  wing  tip  and  compared  to  an 
equivalent  expression  predicted  by  St.  Venant's  theory,  the  fol¬ 
lowing  expression  is  obtained. 
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where  ao  (1)  is  the  wing  tip  twist  given, by  equation  15  while 

aQ  (1)  is  the  wing  tip  twist  given  by  the  St.  Venant  torsion 

theory.  A  plot  of  equation  10  is  shown  in  Fiqure  1. 

5.5.7 

b.  Steady  state  concentrated  tip  twist  loads 

If  the  wing  is  under  the  influence  of  a  concentrated 
twisting  moment,  Fq  at  the  tip  as  a  result  of  a  steady 
state  coupled  bending  torsion  displacement,  the  exact 
closed  form  twist  distribution  that  satisfies  these 
equations  of  motion  and  their  associated  boundary 
conditions  is  given  by 
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41 
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When  the  twist  distribution  given  by  equation  11  is  evaluated  at 
the  wing  tip  and  compared  to  its  counterpart  predicted  by  the 
St.  Venant' s  torsion  theory  the  following  expression  is  obtained. 


ao(1) 

ao(1)  St.  V 


=  1- 


tanh4 A 
_ < 

4A 


(12) 


It  should  be  noted  that  the  ratio  given  by  equation  12  was  plot¬ 
ted  for  the  real  values  of  Ac  in  references  5  and  6  and  was  shown 
to  represent  conditions  where  any  errors  resulting  from  using  St. 
Venant 's  torsion  theory  are  conservative  (over-design  rather  than 
under-design) .  In  this  analysis  equation  12  is  examined  when  Xc 
is  imaginary,  which  is  possible  if  L2  is  very  large.  Under 
such  circumstances,  equation  12  becomes 


ao(D 


a(l) 
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St.  V 
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Figure  2  depicts  the  conditions  given  by  equation  13.  It  is 

therefore  seen  from  the  figure  that  there  are  certain  ranges  of 

Ac  for  which  nonconservative  errors  are  possible  by  using  the 

St.  Venant's  twist  theory. 

5.5.8 

Results  and  Conclusions 

The  results  are  shown  in  Figures  1  and  2.  Figure  1  shows  a  com¬ 
parison  of  the  static  wing  tip  twist  obtained  in  the  present 
study  and  that  obtained  via  St.  Venant's  twist  theory  in  the  pres 
ence  of  statically  distributed  forces  and  low  to  moderate 
coupling.  Figures  2  shows  the  trend  for  concentrated  forces  and 
substantial  coupling.  In  Figure  1  it  is  seen 

that  the  presence  of  coupling  makes  the  errors  of  St.  Venant's 
theory  worse.  This  seems  to  suggest  that  the  more  sophisticated 
theory  is  more  important  for  wings  with  coupling  (e.g.,  wings 
aeroelastically- tailored  using  elastic  cross-coupling) . 

ao{1) 

Figure  2  also  shows  that  nonconservative  errors  (| — ppp  |>1)  are 

°o 

are  possible.  St.V 

Using  Figures  1  and  2,  the  following  conclusions  can  be  summar¬ 
ized:  (i)  ignoring  warping  arbitrarily  using  St.  Venant's  theory 

could  result  in  very  significant  errors  (as  high  as  over  80%  er¬ 
rors)  in  analytical  results  for  composite  aircraft  wings,  (ii) 
warping  is  more  important  (St.  Venant's  theory  is  less  accurate) 

for  wings  with  coupling,  (iii)  St.  Venant's  theory  (which  has  al- 

(12) 

ways  been  shown  to  be  conservative  '  ),  can  be  non-conservative 

or  St.  Venant's  approximation  can  lead  to  an  unsafe  design  error 
(under  design  rather  than  over  design  from  a  stability  point  of 
view)  . 
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Figures 


Figure  1:  Wing  Tip  Twist  Ratios  for  Simple  and  More  Involved 

Theory  (distributed  load  and  low  to  moderate  coupling) 

Figure  2:  Wing  Tip  Twist  Ratios  Comparing  Simple  and  More  In¬ 
volved  Theory  (with  substantial  coupling) 
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6.0  STATUS  OF  PUBLICATIONS 


The  significant  results  obtained  thus  far  in  this  research  oro- 
aram  are  being  compiled  and  written  up  for  publication  in  suit¬ 
able  technical  journals.  AIAA  Paper  No.  86-1006,  a  paper  based 
on  the  first  phase  of  this  research  has  been  prepared  and  was 
oresented  at  the  AI.AA/A3ME/ASCE/AHS  27th  Structures,  Structural 
Dynamics  and  Materials  Conference  held  in  San  Antonio,  Texas, 

May  19-21,  1986.  A  paper  "Exact  Solutions  to  Aeroelastic  Oscil¬ 
lations  of  Composite  Aircraft  Wings  with  Warping  Constraint  and 
Elastic  Coupling"  has  been  submitted  to  the  AIAA  Journal. 

Another  paper  entitled  "Some  Implications  of  the  Warping  Re¬ 
straint  on  The  Behavior  of  Composite  Anisotropic  Beams"  whose 
content  forms  the  basis  for  this  report  has  been  accepted  for 
publication  by  the  Journal  of  Aircraft.  In  addition,  a  paper 
entitled  "Exact  Closed  Form  Solution  for  Nonlinear  Unsteady 
Aerodynamics",  which  provides  the  basis  for  the  proposed  next 
phase  of  the  program,  has  been  accepted  for  publication  by  the 
AIAA  Journal.  Furthermore  the  results  of  the  divergence  study 
are  in  the  process  of  being  written  up  for  submission  to  be  con 
sidered  for  publication  in  a  technical  journal. 
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